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1 . EINLEITUNG 



In der Geodaes i e werden immer haeuf i ger Verfahren und 
Instrument© benutzt, die es erlauben, die Beobachtungsdaten 
ohne irgendwelche Vor ver ar be i tung direkt in den Rechenanlagen 
we i t er zuver arbe i ten . Ein Problem sind auch die grossen 
Datenmengen., die in der Geodaes i e anf a I I en . Dadurch wird es 
immer schwieriger, grobe Fehler in Beobachtungen zu erkennen 
und zu el iminieren. Wie auch al I gerne i n bekannt ist, werden 
durch grobe Fehler (Ausreisser) die Ergebnisse bei der heute 
meist verwendeten Ausgleichung unter Minimierung der 
Quadratsumme der Verbesserungen sehr oft stark verfaelscht und 
verzerrt. Daher sind in den letzten Jahren Verfahren zur 
Aufdeckung grober Fehler entwickelt worden. Die bisher 
verwendeten Methoden der Ausreissersuche nehmen Hypo thesen t es t s 
aus der Statistik zu Hi I f e Cz.B. "data-snoop ing" ) . Aber nicht 
nur das Erkennen grober Fehler in den Beobachtungen ist 
wuenschenswer t , sondern auch das Auffinden grober Fehler in den 
Naeher ungswer t en fuer die unbekannten Parameter (wie z.B. 
Koordinaten usw.), Wenn grobe Beobachtungsfehler und grobe 
Fehler in den Naoher ungswer t en gleichzeitig auftreten, so wird 
eine richtige Diagnose bei Verwendung der Methode der kleinsten 
Quadrate sehr schwierig. 

In der vor I i egenden Arbe i t w i rd e i ne andere Form der 
Ausgleichung behandelt und zwar die Ausgleichung durch 
Minimierung der Absolutsumme der Verbesserungen. Fuer zwei 
geodaetische Anwendungen, 

1 .) Ausgleichung von Lagenetzen 



2.5 He I mer t tr ansfor ma t i on 
werden theoretische und praktische Untersuchungen 
dur chgef uehr t . Es so I I die Frage beantwortet werden, wie sich 
die Abso I u t summenm i n i m i er ung fuer die Suche nach groben Fehlern 
in den Beobachtungen und in den unbekannten Parametern eignet. 



Das Manuskript wurde mittels eines Tex t ed i t or pr ogr ammes 
erstel It und geplottet. Dadurch sind gewisse Abweichungen von 
der deutschen Rechtschreibung (Umlaute, scharfes s) bedingt. 



2. VERSCHIEDENE NORMEN UND IHRE SPEZIELLEN EIGENSCHAFTEN 



Die p-Norm eines n — d i mens i ona I en Vektors v ist in 
f o I gender We i se def i n i er t : 



i v »p = A/ Li v «r 



i= 1 

So I I die Norm unter gewissen Nebenbedingungen zu einem 
Minimum werden, und setzt man fuer p verschiedene Werte ein, so 
koennen daraus un terschiedl iche Mini mumspr inzipien abgeleitet 
werden . 

I . Fuer p = 2 : 

/ i 

2 n 



|v| 2 =A/^ v i = Min 

Diese Form wird 2-Norm genannt und ist gleichbedeutend mit der 
Gauss ' sehen Mini mumsbed i ngung . 

n 

v = Min 



i = 1 



2 . Fuer p = CO 



p' n 



i v| L 

CO p -* oo 



r lim A/ Hl V i| P = Min 



i =1 



Diese Norm ist als Tschebyschef f — Nor m bekannt. Ihre 
Minimierung ist gleichbedeutend dem MINMAX Prinzip 

|| v || = Max |v;| = Min 



3. F 



uer p = 



|v| 1 = £jv.| r Min 

i = 1 



Bei dieser Form wird die Absolutsumme der Elemente des Vektors 
v zu einem Minimum gemacht. 

Sei eine Messreihe von direkten Beobachtungen, nach 
aufsteigender Groesse sor t i er t , gegeben, so zeigen sich bei 
Minimierung der Verbesserungen nach den verschiedenen Normen 
folgende signifikante Eigenschaften. 

Messreihe C 2, 2 , 2 , 2 , 2 , 1 5, 20 ) 



2-Norm : 
ausgeg I i chener Wert — 6.4 (ar i thmet i sches Mittel) 
Verbesserungsvektor CA. 4, 4 . 4 , 4.4, 4.4, 4.4, -8.6, -13.6) 

oo-Norm : 
ausgeglichener Wert = II CMittel aus max. und min. Messwert) 
Verbesserungsvektor (9, 9 , 9, 9, 9, —4, —9) 

Bei der Tschebyschef f -Norm wird nach der Forderung — minimiere 
die maximale Verbesserung — ausgeg I i chen . 

1 -Norm : 
ausgeg I i chener Wert = 2 (Median) 
Ver besser ungs vektor (0,0,0,0, 0,-13,-18) 

Die Minimierung der I -Norm f uehr t auf den sogenannten 
Medianwert. In einer der Groesse nach sortierten Zahlenreihe 
ist es immer der in der. Mitte gelegene Wert. Bei einer geraden 



Anzahl von Elementen eignen sich zwei Werte fuer den Median und 
zwar Element n/2 und Element (n/2)+1 . Auch jede 
dazw i sehen I i egende Zahl kann als Median genommen werden. 



Vergleicht man die nach den verschiedenen Prinzipien 
ausgeg I ichene, mi t zwei groben Fehlern behaftete Messreihe., so 
koennen folgende Seh I uesse gezogen werden: 

Der ausgeg I i chene Wert • wird bei den Verfahren der 2-Norm 
und OO — Norm stark verzerrt. Ebenso erhaelt man durch die 
Verbesserungen keine (OO-Norm) oder nur schlechte (2-Norm) 
Information ueber die groben Fehler in den Beobachtungen. Ein 
befriedigendes Ergebnis erreicht man nur durch die Ausgleichung 
nach der I -Norm . Sie ist fuer eine direkte Beobachtungsreihe 
am besten dazu geeignete grobe Messfehler aufzudecken und auch 
einen vernuenf t i gen ausgeg I i ebenen Wert zu erzielen. Der Grund 
ist in der Robustheit des Med i ans zu finden. Wuerde man zum 
Beispiel in der Messreihe den 7-ten Beobachtungswert (20) 
bei i eb i g vergroessern, so aender t sich der Medianwert dadurch 
nicht. Das Gleiche gi It fuer eine Beobachtung, die kleiner als 
der Median selbst ist. Eine bei i eb i ge Veraenderung nach unten 
bewirkt keine Veraenderung im Medianwert. Man kann daher jede 
Beobachtung veraendern Cunter der Bedingung, dass dabei das 
Vorzeichen der Verbesserung vor und nach der Veraenderung nicht 
gewechselt wird), der Medianwert bleibt fest. In den 
Verbesserungen wird diese Veraenderung durch eine 
Vergr oesser ung des Absolutwertes der Verbesserung angezeigt. 
Es ist leicht einzusehen, dass in diesem Sinne auch mehrere 
oder al le Beobachtungen Causser jener, die dem Medianwert 
entspricht) veraendert werden duerfen. 



FRIEDRICH, K. C 1 937) hat schon die Minimierung von 
Absolutsummen bei der Suche nach optimalen Gewichten in 
Bas i sver gr oesser ungsne t zen verwendet. Die Ausgleichung nach 
der OO-Norm fuer eine He I mer t tr ansf or ma t i on wurde von MEISSL, 
P. C1968) behandelt. Ebenso haben REINHART, E. C 19753 und auch 
HEINDL, G. und E. REINHART CI976) die Ausgleichung minimaler 
Maximalfehler untersucht. BARRODALE, I. C1968) hatte 
numerische Daten analysiert und wies in seiner Arbeit auf die 
Vortei le der I -Norm gegenueber der 2-Norm hin. Die ersten 
Anwender der I -Norm auf geophys i ka I i schem Gebiet waren 
CLAERBOUT y J.F. und F. MUIR C1973). 

Es stel It sich nun die Frage, ob diese Medianeigenschaften 
auch in einem al (gemeineren Model I bestehen bleiben und wie man 
die konventionel I en Model le in der Ausgleichsrechnung unter der 
Forderung - minimiere die Absolutsumme der Verbesserungen - 
berechnen kann. Die Loesung auf diese Fragen bietet sich in 
der LINEAREN PROGRAMMIERUNG C lineare Optimierung) an. Daher 

werden zuerst im naechsten Kapitel einige grundsaetzliche 
Theoreme der I inearen Programmierung besprochen. 



GRUNDLAGEN DER LINEAREN PROGRAMMIERUNG 



3.! Formul i erung des mathematischen Model \s- 

Bei der Formul i erung des mathematischen Problems wird 
zwischen der Standardform und der Norma I f orm unterschieden. 

3.1.1 Standardform eines linearen Programmes: 

a „y, + a i2 y 2 + • - • + a nyi - b i 

C3. I ) 

a mi y, + a m2 y 2 + ■ • ■ + Q ml y. - b m 

c, y, + c 2 y 2 + • • + c ( y ( = Max zieifunktion 

y, > 0, y 2 ^ 0, ■ ■ J^O 

Sollte eine der Ungleichungen mit der Relation >= auftreten, so 
kann durch Multiplizier en mit — 1 die Relation umgekehrt werden. 
Weiters ist leicht einzusehen, dass durch Mu I t i p I izieren der 
Zielfunktion mit —1 die Forderung nach Minimieren der 
Zielfunktion erreicht wird. Die zweite wichtige Formul ierung 
ist : 

3.1.2 Normalform eines linearen Programmes: 

Die Normaiform eines I inearen Programmes wird durch Umformen 
der Ungleichungen in Gleichungen bewer ks teil igt. Dies wird 
erreicht durch Einfuehren von zusae tz I i chen ni chtnegat iven 
Variablen. Diese Variablen werden Seh I upf var i ab I e genannt. 
D i e Ung I e i chungen aus 3.1.1 erha I ten dann f o I gende Form : 
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a„ y, + • ■ • + a„ y, + y, 



= b, 



amiyi 
c, y, 



+ a m , y, + 
+ c, y, 



= Max 



C3.2) 



yi -°, • • , y< -0, y, * 0, • • • , y m >o 



D i e Var i ab I en y , . . , y s i nd d i e schon vorher erwaehnlen 

1 m 

Seh I upfvar i ab I en. C3 . 2) stellt bereits eine spezielle 
Auspraegung der Normal form dar. In weiterer Folge wird nur 
mehr die al I gerne i ne Form der Norma I f orm betrachtet. Es werden 



alle Var i ab I en y , . . . , y ; y 
zusammengef asst . 



' Y« 



zu einem Vektor x 



a,, x, + 



+ a !n x n 



b, 



m < n 



a m) x. + 



C ! *, + 



+ ^mn X n - D m 



• +■ c n x n = Max 



C3.3) 



x, >0, • -, x n >0, b, >0 , • • -, b m >0 



3.2 Die Simplexmethode 



Ist in den Gleichungen C3.3) kein Widerspruch vorhanden, 
so gilt wegen m<n eine n — m par ame tr i ge 
Loesungsmann i gf a I t i gke i t . Man koennte zum Beispiel nach den 



ersten Unbekannten 



Ci = i, . . 



au 



f I oesen Cunter der 



Annahme,, dass die entsprechende Un t er de t er m i nan t e von Nu I I 
verschieden ist) und al I en restl i chen Unbekannten 
x m*1 ' x m*2'''-' x n c ' en Wert null zuwei sen . Wir nehmen dabei an, 
dass x >= Ci = I , . . . , nO ist. 

D i e Unbekannten x, ,x v . . w x m >= werden dann 
Basisvar iable y die restlichen x^ , , . . . , x „ N i ch tbas i svar i ab I e 
genannt. Es ist leicht einzusehen,, dass natuerl ich auch eine 
andere Kombination von Basis— und N i ch tbas i svar i ab I en benutzt 
werden kann. Voraussetzung ist jedoch, dass jedesmal das 
Gleichungssystem nach den Basisvariablen auf loesbar sein muss 
und dass sich dabei bei verschwindenden N i ch t bas i svar i ab I en 
nichtnogat i ve Werte fuer die Basisvariablen ergeben. 

Die grundlegende Idee des Simplexalgorithmus ist es nun, 
durch sukzessives Austauschen der Bas i svar i ab I en mit den 
N i ch t bas i s var i ab I en die Zielfunktion zu maximieren. Das 
Ausgangsstadium fuer jeden Simplexschritt ist die kanonische 
Form. Diese Form wird durch Ausdruecken der Basisvariablen 
durch die N i ch t bas i svar i ab I en in den Zwangsbedingungen und in 
der Zielfunktion erreicht. 

Angenommen, es ist eine zu I aess i ge Basisloesung von C3.3) 
vorhanden, das he i sst alle Bas i svar i ab I en s i nd >— . Di ese 
Basisloesung bei nha Ite zum Beispiel x. >— 0,..., x m >~ als 
Bas i svar i ab le und x mvl , . . . , x n als Ni ch tbas i svar i.ab I e . Dr ueck t 
man die Basisvariablen durch die N i ch tbas i svar i ab I en aus, so 
erhaelt man die kanonische Form. 
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+ a, 



X m + 1 + 



+ a 



in 



b, 



*m + a m.m»1 *m.| + 



BASISVARIABLE 



+ a „ x„ 



NICHTBASISVARIABLE 



= b 



C3.4) 



Z + C m*1 X m.1 + 



n n 



Die Koeffizienten c" werden relative Kostenfaktoren 
genannt, z ist der Wert der Zielfunktion fuer die 
augenblickliche Loesung. Die Optimal itaet einer Basisloesung 
wird durch die relativen Kostenfaktoren c . Cj = m+ 1 , . . . , rO 
angeze igt. 

Sa tz • 

Eine zulassige Basisloesung ist maximal Cop t. i ma I 3 , falls alle 

c >~ sind. Bzw . y wenn die Zielfunktion zu minimieren ist, 

wenn alle c. <= s i nd . 
J 

Bewe i s : 

Der Beweis ist hier nur fuer ein Maximieren der Z iel funkt ion 
gefuehr t. Sei x>x 2 >.>>< n einebeliebi ge andere zulaessige 
Loesung Cm i t x >= 0, . . . y x >= 3 0) und se i en i n der kanon i sehen 
Form C3.4) alle c. >= Cj = m+] , . . . , n) , so ergibt sich ein 
neuer Wert der Zielfunktion z Q fuer die Loesung x , x ? , . . ,x 



^o ~ ^o ^m«-l *m*1 



- C X < 7 

^ n A n — *-o 



Wegen c , >= 0, . . . , c >= 
^ m + 1 n 



m* 1 



>= 0. 



x >= wird das 

n 



vorhergehende z o auf jeden Fa I I > — z Q sein. Das heisst y die zv 
z o gehoerende Basisloesung x. , . . . , x n ist optimal. 



Waere nun z.B. der Z i e I f unk t i onskoef f i z i en t c < 0, so kann 



durch Wah i vo 



n x , 



> die Zielfunktion z vergroesser t 



werden. Damit wird aber x Basisvariable. Wie gross darf nun 
gewaeh I t werden? Nach der N i ch t nega t i v i t ae t sbed i ngung darf 
keine Variable negativ werden. 
Aus (3.4) ergibt sich 



x, = b, - a, s -G 



x m = b m - a ms Ö 



Damit keine der Variablen negativ wird, muss 



= Min -J*- ( 1 = 1,2, • • .,m ) , a is > 



gewaeh I t werden. Dieser Vorgang wird auch Pivotsuche genannt. 
Sei nun r jener Zeilenindex, fuer den das Minimum eintritt^ so 
i s t 



= 



cu 



und die Unbekannten lauten: 
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b. - a, ■— - 

1 IS Q 



x r = b f - a rs -=!- 







m m ms q 



b; 



x s = 



Wie aus den neuen Unbekannten zu ersehen i st, wurde x zu nul I 
gemacht und x s > 0. Diese Operation wird Basiswechsel genannt. 
Es wurde die Spalte der Unbekannten x in die Basis genommen 
und die Spalte der Unbekannten x aus der Basis geholt. Dieser 
Basiswechsel wird sukzessive so lange dur chgef uehr t , bis alle 
c > = sind, das he issl, das Optimum erreicht ist. Minimiert 
man die Z i e I f unk t i on, dann ist die Loesung optimal, wenn alle 
c <- s i nd . 

Aus dem bisher Erwaehnten ist leicht zu erkennen, dass 

jede Bas i s I oesung auch durch Auf I oesen eines m mal m 

Gleichungssystems zu erhalten ist. Die m Spalten dieses 

Systems sind ident mit den jeweili gen Spalten der 
Bas i svar i ab I en . 



Die Theorie der I i nearen Programmierung I i ef er t auch den 
Beweis fuer die Existenz einer optimalen Loesung nach endlich 
vielen Simplexschritten. Auf die Beweise dazu, sowie das 
Auftreten von Zyklen und entarteten Loesungen sei auf DANTZIG 
G.B. C1963) und COLLATZ L. C1971D verwiesen. 
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Ein wichtiger Abschnitt der linearen Programmierung 
behandelt das Dual i t ae t spr i nz i p . 

3.3 Dualitaet bei Normalform 



Wird die Norma I f orm C3.3) 



n 

n 

) c; Xj = z = Max 



I 



C3.5) 



M 



c T x = z- Max 

als pr i ma I es Problem betrachte t y so kann ein duales Problem 
gegenueber ges t e I I t werden und zwar in der Form 



Ml 



i l 



TT. > c 



1=1 

m 



^b, IT; = u = Min 



I 



C3.63 



1 = 1 



A T I>c 

b U = u = Min 



TT; 



'orze i chenunbeschr aenk t 



i = 1 , ■ . . y m ; j ■ I , . . . , n 



D i e vorze i chenunbeschraenk ten Unbekann ten U". Ci = I , . . . , mO 
des Dualen Problems werden auch Multiplikatoren genann t . Im 
folgenden werden zwei wichtige Saetze der Dual i taet angefuehrt . 
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Die Beweise dazu werden hier nicht gebracht, sie sind aber in 
den meisten Buechern ueber lineare Programmierung behandelt 
(siehe z.B. DANTZIG G.B. C1963) oder auch COLLATZ L. C197135. 

Satz: 

Fa I Is CD und CID zu I aess i ge Loesungen haben, dann gi It fuer 

die zugehoerigen Werte der Z i e I funk t i or 

2 <= U 

MAXCz) <= MINCu) 



Sa tz : 

Falls CD eine Op timalloesung hat, so besitzt auch CID eine 

so I che und es gilt 



z = u 



bzw . 

MAXCz) = MINCu) 



Ausserdem existieren noch verschiedene Regeln, die die 
Zusammenhaenge zwischen pr imalem und dualem Problem 
beschreiben. Sie werden Kuhn-Tucker Regeln genannt. In 
folgender Tabelle sind diese Regeln dargestellt. Man liest das 
Primale zeilenweise und das Duale spaltenweise. CD soll 
andeuten, dass die dazugehoer i ge Variable keiner 
Vor ze i chenbeschr aenkung unter I legt. 





PRIMAL 






VARIABLE 


X, > • 


• • x.(t) . • 


■ V^ 


RELATION 


KONSTANTE 


n,>o 


a„ • 


• . Q,j . ■ 


• a, n 


< 


b, 




. 


'. 


• 


* 




. 




• 




• 


• 




• 


_J 














< 


TTi(i) 


a-,i ' 


• • ^ ■ • 


* a. 

in 


— 


b-, 


o 


• 


. 


• 


• 




• 




• 


. 


• 


• 




• 




• 


• 


• 


• 




• 




iL*o 


Q m, • 


- ' Q mj ' 


. . a 

mn 


< 


b 

m 


RELATION 


> 


= 


> 




i 

.E 












KONSTANTE 


c, • 


. . Cj .. ■ 


• c n 


— Max 


2 



Tabel le 3. 1 
Kuhn-Tucker Rege I n 
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Es ergeben sich daher folgende anzuwendende Regeln 



PRIMAL 


DUAL 






Koeff izi entenmatr ix A 
Zioifunkt ion MAX 
constante G i ieder 


transponierte von A 
constante Gl i eder 
Zie 1 funkt ion MIN 


Vor j ab 1 o : 
x. >« 

x- keine Vorzeichen— 

beschraenkung 


Rq 1 a t i on : 

Cj—te) Ungleichung > — 

Cj— te) Gleichung = 


Re 1 a t i on * 

Ci— l©) Ungleichung < = 

vi— te) Gleichung - 


Var i ab 1 e * 

Ci-ter) Multiplikator >= 
C i - t er ) Mu Itiplikator keine 
Vorze i chenbeschr aenkung 



An Hand dieser Regeln und der vorangegangenen zwei Saetze ist 
es leicht einzusehen^ dass das Duale von CID wiederum das 
pr imale Problem CID ergibt. 



Im folgenden ist noch der Satz vom komplementaeren Schlupf 
zu behandeln. Seien das Primale und Duale als Systeme von 
Ungleichungen CS t andar dform D vorgegeben und durch Einfuehren 
von Seh I upf var i ab I en in Gleichungen umgewandelt, dann lautet 



das Pr i ma I e : 
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a„x, ♦ 

a 21 x,* 



QmlX, + 



C, X, * 



•1 J 



• + a 2 n X n 



+ a mn x n 



+ x 



n * 2 



+ c n x n 



♦ x n . 



= b. 



= b, 



"rul J 



x n >0 , b n • • ■ b m >0 
. x £ ■ • ■ Schlupfvariable 



Ma x. 



und das Dua I e 



a H TT, * 
a 12 TT, ♦ 



+ a m .TL - TT m , 



•■m.'m 



m ♦ 1 



TT 



m-2 



= c, 

= C, 



b, TT, * 



+ a TT 

mn m 



b m n 



m u m 



TT 



= Min. 



TT,, - TT m *0 

TT •TT ^0 



Schlupfvariable 



m i i • i ■ . i i -7i i C i = 1 , . . . j n) des Dua I en m i t der 

Mu I t i p I i z I er t man jede Ze i le ' 

korrespondierenden Unbekannten x . <« = 1 n> des Pr i ma I en, 

summiert diese Zei len und subtrahiert sie von der Ziel funklion 
des Dualen, sa erhae I t man bei einer optimalen Laesung 



n 

(bi .^a,,x,)Tr, +(b 2 - IIa 2l x i )TT 2 ♦ • ■ • ♦ l b m - Y1^^)K 



**lK.l + X 2 1T m.2 + 



= 1 

+ x TT 

n m ♦ n 



= 



Auf der rechten Seite dieser Gleichung steht die Differenz der 



8 



Z i e I f unk t i onswer t e der pr imal en und dualen Loesung . Ist die 
LoGsung optimal, so sind beide Ziel funk t ionswertö g I e i ch / daher 
die Null auf der rechten Seite. 
Die Kl ammer ausdr uecke 



ii 

(b j ■ IZ a ii x > } j= ] >- ■ -i m 



sind aber die Seh I upf var i ab I en x, Ck 



j n+m) des 



Pr imal en. Jedes Glied auf der linken Seite ist ni ch tnega t i v 
(wegen x. >= und TT. > = Ci = I , . . . , n+m) ) , das heisst jedes 
Gl i ed muss einzeln verschwinden. Daher muessen folgende 
Beziehungen gelten: 



1 1 

(bj - ^Zdj.x, )TTj = j= 1 r • -m 

■■ k =1 • • • n 



i = i 



k m ♦ k 



Daraus kann nun geschlossen werden: 

Befindet sich eine Variable des dualen Problems in der 
Basis Cund zwar handelt es sich um die Basis einer optimalen 
Loesung), so ist der Schlupf der entsprechenden Ungleichung des 
Pr i ma I en gleich nu I I Cund umgekehrt,, da ja das Duale wieder als 
Primales und das Primale wieder als Duales aufgefasst werden 
kann). Weiters muss fuer jene Beziehungen des Pr i ma I en 
(Dual en) mit der Relation < (>) die korrespondierende 
Unbekannte des dualen (pr i ma len) Problems nu I I sein. 



Es gibt nun verschiedene S i mp I exa I gor i t hmen , die die 
Multiplikatoren mitliefern. Einer davon ist die revidierte 
S i mp I exme t hode , die im folgenden Kapitel kurz beschrieben wird. 
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3.4 Revidierte Simplexmethode 

Man betrachtet wieder die Normalform C3.3) und nimmt an, 
dass eine zu I aess i ge Bas i s I oesung vorhanden sei . Die 
vorhandene Bas i s I oesung ist z.B. 



x , =q, 



X 

m 


zz 


q m 


X 


zz 


0" 


rrw 






X 

n 


zz 


0_ 



-Basis 



Nichtba&is 



Die Spalten der Koeffizientenmatrix A der Normalform - geb i Idet 
aus den Basisvariablen — koennen zu einer Matrix B 
zusammengef ass t werden. Mit der vorhin als bekannt 
vorausgesetzten Bas i s I oesung sieht dies f o I gender massen aus: 



ia n x, + 



D 



a m i x . + 



c, x, + 



+ a, x i+ a, , x , + 



+ Q mm X m 



c^x 



*" a m.m*l X m*1 * 



+ C 



m*1 A nru1 



BASIS 



+ Q 1n X o 

i n n 



*a m n X n 



+ c n x n 



NICHTBASIS 



= b. 



= z = Max 



C3.75 



Die In der Normal form C3.7) eingezeichnete Matrix B und der aus 

den Kos t enf ak toren Cj Ci = 1 , . . . y m) gebildete Vektor c liefern 

die Multiplikatoren TT- Ci = ! , . . . , m) . Sie sind die Loesungen 

der G I e i chungen 
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B T I = c 

Mit Hilfe der Multiplikatoren TT erhaelt man auch die 
relativen Kostenfaktoren. Werden saemtl i che Gleichungen der 
Normalform C3 . 35 mit don Multiplikatoren [T Ci = 1 y . . . y rrO 
mu I t i p i i z i er t , auf summ i er t und von der Ziel funk t i on 

subtrahier t y dann ergibt sich: 



■ ■ ■ in 



i n 

^„-^TtT; a in )x n 



i = l 



Ul 



u rru 1 



n 

= d ♦ Yji x j 



j = m -1 



In dieser Form ist d der momentane Wert der Zielfunktion fuer 



d i e 



zu I aess i ge 



B 



as i s i oesung 



•x, = q, * 



* x, 



Qm- 



x . — ..... x = 5. und die c entsprechen wieder den relativen 
Kostenfaktoren wie in Formel C3.45. AI lerdings sind sie hier 
mit dem entgegengesetzten Vorzeichen behaftet. Sie werden aber 
in der Folge immer auf diese Weise definiert sein. Das heisst, 
das Optimum ist errei cht , wenn alle c • Cj = m+1 , . . . , n5 kleiner 
gleich nu I I sind. Es wird nun gleich wie bei der 
Simplexmethode verfahren. Sei nun irgendein c« > 0, so ist 
eine Verbesserung des Wertes der Zielfunktion z moeg 1 ich. Ist 
c s > y so kann wiederum 



x s =0 > 



gesetzt werden. Wie gross darf man waehlen? Aus den 
Gleichungen C3.7) wird erhalten 



/ | 
X, 



B 



'»,» 



1 b m 



a, 5 



a ms / 



•0 



B" 



Multipl iziert man von rechts mit B " ' (B ~ muss ja fuer j ed< 



zu I aess » 



i ge Basisloesung vorhanden se in), so gilt 



\ X r 



1 qj 



/ t,\ 



i L/ 







w 



iederum erhaelt man 9 mit 



= Min 



tj>0 



ti 



tr 



Der neue Wert der Zielfunktion ergibt sich dann mit 



z = d + c 5 0=d>d 



Dies entspricht wieder einem Basiswechsel . Es wurde die 

genommen,, bzw. x r ist 
N i ch t bas i svar i ab I e geworden. Die neuen Unbekannten lauten dann 



Unbekannte x 5 in die Basis 
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x, = q, - t, = q, 



x r q - 1 ,0 = 8 , 

m-1 n m-t m-1 ~m-1 



x s =0 



-BASIS 



J 



x m =0 



x r =0 



x n =0 



•NICHTBASiS 



Bei diesem Basiswechsel wurde aber nur eine Spalte der 
Basismatrix B geaendert. Verfqhren, welche die neue Inverso 
einer Matrix berechnen, an der nur eine Spalte ausgetauscht 
wurde, werden in der schon vorhin erwaehnten Literatur CDANTZIG 
G.B. CI9635, COLLATZ L. C 1 97 1) ) beschrieben. Ebenso kann dort 
der genaue Algorithmus fuer die revidierte Simplexmethode 
nachgelesen werden. 

Das Optimum ist dann erreicht, wenn al le relativen 
Kostenfaktoren c der Basisspalten nu I I sind, jene der 
N i ch t bas i sspa I ten <- sind. Der Wert der Zielfunktion ist 



m in 



i = i 



i =1 



e i n 



Man sieht, dass im Prinzip bei jedem Basiswechsel 

G I e i chungssys t em Bx = b mit B (m«m) zu I oesen ist. Die 

Spalten von B werden aus den Spalten von A gebildet, die zu den 
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jeweiligen Basisvariablen gehoeren. Die revidierte 
Simplexmethode hat den Vor te i \, dass man immer nur die 
Basismatrix B mit einer neuen ausgetauschten Spalte berechnen 
muss . Ausserdem kann man immer auf die urspruengl ich 
vorgegebenen Werte der A Matrix, der rechten Seite b und der 
Z i e I f unk t i onskoef f i z i en ten zur ueckgr e i f en . Die Simplexmethode 
I iefert ausser der Loesung x auch noch die Mu Itiplikatoren. 

Ein wichtiger Satz der Theorie der I inearen Optimierung 
ist noch zu erwaehnen . Er wird fuer spaetere Betrachtungen 
sehr wertvol I sein. 

Satz : 

Werden bei einer optimalen Bas i s I oesung Ausgangsparameter Cwie 
zum Beispiel die rechte Seite u.s.w. ) geaendert, dann bleibt 
die Loesung optimal, solange sie zu I aess i g bleibt, und solange 
die relativen Kostenfaktoren ihr Vorzeichen nicht aendern. 
CAendert man nur die rechte Seite, so bleiben die 
Kostenfaktoren gleich. Es ist dann nur die Zu I aess i gke i t zu 
f or der n . ) 

Fuer ein genaueres Studium der Simplexmethode und der 
revidierten Simplexmethode ist wieder auf die Buecher von 
DANTZIG G.B. CI963) und COLLATZ L. (197!) verwiesen. Diese 
komprimierte Zusammenfassung diente nur der Erklaerung einiger 
grundlegende Begriffe, Bezeichnungen und wichtige Saetze, die 
spaeter noch Verwendung finden. 
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Weiters Ist hier noch zu bemerken, dass das Auffinden 
einer zu I aess i gen Anfangsbas i s I oesung hier nicht beschrieben 
wird. Verfahren dazu sind wieder aus der vorhin erwaehnten 
Literatur zu entnehmen. Es soll aber bereits hier darauf 
hingewiesen werden, dass eine zulaessige Anfangs I oesung in den 
h i er vor I i egenden Prob I emen ohne Schw i er i gke i t gefunden werden 
kann . 



4. DAS MODELL DES VERMITTELNDEN AUSGLEICHES UNTER 
MINIMIERUNG DER ABSOLUTSUMME DER VERBESSERUNGEN 



4 . 1 Formu I i erung als I i neares Programm 



Es stellt sich nun die Frage, wie man beim Mode II des 
vermittelnden Ausgleiches die Absolutsumme der Verbesserungen 
minimieren kann. Wie schon vorher erwaehnt, I i eg t die Loesung 
in der I i nearen Optimierung. 

Gegeben ist der vermittelnde Ausgl e ichsansatz : 



bzv 



^x - v - I C4. 1 ) 



A Cm- rO I " Cn - 1 ) 

v Cn * 1 ) x Cm»l) 

unter der Forderung — minimiere die Absolutsumme der 
Verbesserungen. Es bedeuten: 

A ... bekannte Matrix der Koeffizienten der 

I i near i s i er t en Verbesserungsgleichungen 
I ... Vektor der Messungen (gemessene Beobachtung 
minus Naeher ungswer t der Beobachtung) 
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v ... Vektor der Verbesserungen 

x ... Vektor der unbekannten Parameter 
x,v, I sind hier vor ze i chenunbeschr aenk t und entsprechen daher 
nicht der N ichtnegat i vi taetsbed i ngung aus der linearen 
Optimierung. Um dieser Forderung gerecht zu werden,, f uehr t man 
neue Unbekannte ein, indem man x als Differenz zweier 
nichtnegativer Werte darstel It. 



Xj = Xj - x } 



V . = v. 

I I 



( j = 1 1 ■ • -, m ) 
( i = 1 , • • •, n ) 



Nega t i ve 



I j . . . , n) ko ennen dur ch Mu Itipl izieren der 



Gleichungmit -1 in positive I ■ um ge wände I t wer den 
G I e i chung C4 . 1 ) erhae I t daher d i e Form 



Ax + - A x* - v* + v" = 1 



Die Ziel funkt ion lautet 



) I v. I = M i n 

i = l 



und ist nicht linear. 

Fuehrt man die schon oben erwaehnten neuen v- ein> so kann die 

nichtl ineare Zielfunktion durch eine I ineare ersetzt werden. 

Die beiden Z i e I f unk t i onen stimmen ueberein, wenn fuer jedes i 

entweder v* «= oder vT = ist. Dies ist bei einer 
i i 

allgemein zulaessigen Loesung nicht garantiert. Weiter unten 
wird gezeigt, dass diese Bedingung bei einer optimalen Loesung 
i mmer e i ngeha I ten w i rd . 
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L v ; +v i 



= Min 



i= i 



Das he i sst, der vermittelnde Ausgleichsansatz mit 
Minimierung der Absolutsumme der Verbesserungen praesentiert 
s i ch als lineares Programm der Gestalt 



Ax + -Ax~ - V* + V~ = I 



C4.2) 






♦ v" = Min 



1 = 1 



In diesem I inearen Optimierungsproblem sind al le geforderten 
Restriktionen wie die N i ch t nega b i v i t ae t sbed i ngung in den x,v 
und I er f ue I I t . Es koennte daher mit ei nem der vorher 
besprochenen Verfahren, wie Simplexmethode oder revidierte 
Simplexmethode geloest werden. 

Es laesst sich leicht beweisen, dass immer nur ein v + oder 



v . 



bzw 



x . 



oder 



Ci - 1 , . . . , n ; j — 1 „ . . . , m) mit pos i t i vem 



» J J 

Wert Caiso > 0) auftreten darf. Wird in der Zielfunktion nur 

ein Paar von Verbesserungen, zum Beispiel von i" „ betrachtet, 

und i s t 



v\= 10 

l 

v: =12 



d.h. v. = -2 



dann ist der Antei l y den beide zum Z i e I f unk t i onswer t I iefern, 
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V m + V 
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Da aber die Zielfunktion zu einem Minimum werden so I I, gibt es 
auf jeden Fa I I eine andere Kombination der beiden v.« , die 
unter Einhaltung der Nebenbedingungen einen kleineren Antei I an 
dem Wert der Zielfunktion I iefert y naeml ich 



v- =0 



»;. = 2 



Dies beweist, dass nicht beide Verbesserungen v* und v~ 

i ' i * 

positiv sein duerfen. Eine zweite Ueberl egung, die auf 
dasselbe Ergebnis fuehrt, ist folgende. Angenommen; es sind in 

einer Bas i s I oesung bei Formel (3.7) ein x + und ein x~ 

J* j * 

vorhanden. Damit sind aber in der Bas ismatrix B zwei linear 

abhaengige Spalten Cdie Spalten der beiden x unterschei den sich 

ja nur dur ch das Vorzeichen) vorhanden, und damit ist das System 

nicht nach den Basisvariablen auf I oesbar . Das Analoge gilt 

fuer ein v* und v" Zusammenfassend heisst di es, dass nur 
i • i • 

jeweils eine der Groessen x* oder x~ bzw. v* oder v~ als 

J J j i i 

Basisvariable auftreten darf. Die andere Groesse muss dann als 
N i ch t bas i svar i ab I e aufscheinen. 



Wuerde nun das I i neare Programm (4 . 2) mit ei nem der b i sher 
besprochenen Algorithmen geloest werden, so muesste man bei der 
Simplexmethode eine C2(n-*-rrOD mal n Matrix abspeichern und auch 
bei der revidierten Simplexmethode haette man mit ei n&r n ma I n 
Matr i x zu oper i ^r^n . Ausserdem wuerde nie die spez i e I I e Form 
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der Zielfunktion Cin den Spalten der x* bzw. x7 ist der 

J J 

Zielfunktionskoeffizient null, in den Spalten der v * bzw. vT 

i t 

ist der Koeffizient - I ) ausgenutzt werden. 

BARRODALE I. und F.D.K. ROBERTS C I 973) haben e i netn 
Algorithmus entwickelt, der eine 1-Norm Approximation eines 
ueberbestimmten linearen Gleichungssystems leistet. Dieser 
Algorithmus bietet Vorteile beim Abspeichern der notwendigen 
Daten und ar beitet schneller als die konventionellen 
S i mp I exmethoden . 



4.2 Fa I I direkter Beobachtungen mit graphischer 
Darstel lung des I inearen Programmes 



Es laesst sich leicht beweisen, dass ein vermittelnder 
Ansatz fuer direkte Beobachtungen (unter Minimierung der 
Absolutsumme der Verbesserungen) und die Loesung ueber die 
I ineare Programmierung auf den Medianwert f uehr t . 

An Hand einer Reihe direkter Beobachtungen sei dies 
illustriert. Die Beobachtungsreihe lautet: 

Cl, 2, 3, 4, 7) 

Die Restriktionen des I inearen Programmes sind 



x - v. = I. C i = 1 , . , . , n) 

i t 
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wobei n die Anzahl der Beobachtungen angibt. 

Legt man ein kartesisches Koordinatensystem zugrunde, mit 
v als Ordinate und x als Absz isse, so koennen die Restriktionen 
als Geraden dargestellt werden. Diese Geraden sind in Fig. 4.1 
mit g . . . w g bezeichnet. Alle dieser Geraden er f ue I I en die 
G I e i chungen 



v *= x 



und haben die Steigung 50 Neugrad und den Abstand I vom 
Ursprung. Vor ze i chenbeschr aenkungen in v, x und I sind nicht 
no t wend i g . 



Bi Idet man die Absolutwerte der v in den Geraden g , . . . g , 
so muessen die Geradente i le der negativen v um die x— Achse 
gespiegelt werden. Es entsteht dadurch ein Knick von 100 bei 



jeder Geraden auf der x~ Achse beim Messwert 



N 



un Koennen 



al le absoluten v Werte addiert werden. Man sieht, dass die 
Neigungen der geradl i n i gen Abschn i tte der Summe umso groesser 
wird, je weiter man sich von der "Mitte* der x. Werte entfernt. 
Daraus folgt aber, dass das Minimum der Summe der absoluten v, 
bei jenem x liegt, welches von allen am "meisten in der Mitte'* 
ist, also fuer das kongrete Beispiel in Fig. 4.1 bei x - 3 . 
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Fig. 4.1 
graphische Darstol 1 ung des I inearen Programmes 
bei direkter Messreihe 
(ungerade Anzahl von Beobachtungen) 



Bei einer geraden Anzahl von Messungen kann der Median 
jeden Wert zwischen den beiden Beobachtungen n/2 und Cn/2 + 1D 
annehmen. P i g . 4.2 veranschaulicht die Situation bei gerader 
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Anzahl von Beobachtungen mit der Messreihe 



<l .2.5.8) 



5^ | v [ = Min für 2 < x < 5 




Fig. 4.2 

graphische Darsiel lung bei gerader Anzahl von 

Beobachtungen 
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5. ALGORITHMUS VON BARRODALE UND ROBERTS ZUR 
ZUR I - NORM APPROXIMATION EINES UEBER - 
BESTIMMTEN LINEAREN GLEICHUNGSSYSTEMS 



Eine I ineore I -Norm Approximation geht von folgender 
Aufgabenstel I ung aus. Gegeben ist ein ueberbes t i mm t es, 
inkonsistenles, I ineares Gleichungssystem 



A x = l 

A ( n*m ) 

x ( m* 1 ) 

1 ( n-1 ) 



C5. I 5 



Die Aufgabe besteht dar in, die beste 1— Norm Approximation zu 
suchen, d i e 



n m 



r ir-, 



1*1 jxl 



J I 



minimiert. Dies ist die gleiche For der ung,, die ein 
vermi ttelndes Ausgl ei chsmodel I unter Minimierung der 
Absolutsumme der Verbesserungen stellt. Bringt man wie ueblich 
an den inkonsistenten Gleichungen C5. 1) Verbesserungen v. an 
Csie entsprechen den Seh lupfvar i ab I en in der linearen 
Opt imierung), so so I I 



n n m 



i*i 



i*l j«l 
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minimiert werden. 

Oder analog der urspruengl i chen Formu I ierung 



A X - V = 



1 1 

Y |Vi| = Min 



1 = 1 



Die durch E i nfuehren von neuen Unbekannten x. = x + — x und 

J J J 

v. - v. — v 7 erhaltene Form (4.23 kann in e i nem Tab I eau 
i i i 



dargestel It werden. 



^\ 


-Ä 


. =D ' 


r— | 



o 



o 



-1 



Max 



Fig. 5.1 
Tab I eau zum Algorithmus von Barrodale und Roberts 



In diesem Tableau wurden etwaige Zeilen mit negati ven I schon 
vorher mit — 1 multipliziert. Da dieses Tabl eau 'dem 
S t andar d t ab I eau der Simplexmethode entspricht, wo die 
Zielfunktion maximiert wird, f uehr t man negative Vorzeichen in 
den Ziel funkt i onskoef f iz i enten g \ n Cdaher die —1 in d^r 
Z i e I f unk t i on ) . Betrachtet man dieses Starttableau etwas 
genauer , so kann man erkennen, dass eine zu I aess i ge 

Anfangsbas i s I oesung sofort vorhanden ist. In dieser 
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v ? Ci - I , . . . , n) in der Basis. -Wird 



Anfangs I oesung sind alle 
die Anfangs I oesung v~ — I in die Zielfunktion eingesetzt, 
erhae I t man die relativen Kostenfaktoren fuer die Anfangsbasis. 
Das Tab I eau in Figur 5. I bekommt dann die Form 



^\ 


■2^ 


a D 


D 



I =1 ul 



"7i fi 

ß-i li A- rni 

i=1 i = l 











-D 



i=i 



Fig. 5.2 
Tab I eau mit Anfangs I oesung und dazugehoer igen 
relativen Kostenfaktoren 



Die relativen Kostenfaktoren bestehen also bei den 
Unbekannten x* , x Csie sind N i chtbas i s) aus der Summe der 
dazugehoer i gen Spalten und aus — 2 bei den 
Nichtbasisverbesserungen. Wie man auch sieht, muss die Summe 
der relativen Kostenfaktoren zweier korrespondierender x* und 
x " null sein, die Summe der relativen Kost enf ak tor en zwe i er 
korrespondierender v * und v * muss -2 ergeben. We i ters kann 
leicht erkannt werden, dass nur die Matrix A, der Vektor \ und 
ein Vektor fuer die relativen Kostenfaktoren, ferner zwei 
Hi I fsvektoren Cfuer Basis— und Nichtbasisi ndizes) zu speichern 
sind, um die vo I le Information ueber das Problem zu haben. Das 
heisst, das Ausgangs t ab l eau kann in der Form, wie in Figur 5.3 
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dargeslel I t .. gebi I de t werden 



INB 




CO 



Fig. 5.3 
Ausgangslabl eau mit Hilfsvektoren 



In Figur 5.3 bedeuten 



9 

A 

L 

INB 

IBASIS 

C 

CO 



Koef f I z i entenmatr i x 

rech t e Se i t e 

Vektor fuer N i chtbas isvar i ab I e 

Vektor fuer Bas i svar i ab I e 

relative Kostenfaktoren 

Wert der Ziel funkt ion fuer die augenb I i ck I i che 

Basi sl oesung Cbe i der Anfangs I oesung idont mit 

der Summe der Gl i eder der rechten Seite) 



Der Algorithmus arbeitet in zwei Stufen. Stufe 1 
beschraenkt die Wahl einer Pivot — Spalte Calso einer 
Unbekannten, die in die Basis kommt) waehrend der ersten m 

Cj = 1 < . . . , fiO . Der 
Vektor, der in die Basis kommt, wird gewaeh I t als jener mit dem 



Iterationen auf die Spalten von x? und x 7 

J J 
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groessten positiven relativen Ziel funkt i onskoef f i z i enten. Der 
Vektor der in die Nichtbasis kommt, wird aus den Vektoren v + 
und vT gewaehl t . Und zwar jener, der die groesste Reduktion 
des Z i e I f unk t i onswer t es bewirkt. Am Ende von Stufe I ist der 
Rang k C<=nO der Matrix A bestimmt. Das momentane 
Simplextableau repraesentiert eine Approximation, die genau k 
Gleichungen \^on CS. 2} konsistent macht. In Stufe 2 wird nur 

mehr getauscht zwischen Basisvektoren v * oder vT mit 

» i 

N i ch t bas i svek t or en v* oder vT. Basisvektoren x* oder x" ist es 

11 J j 

nicht erlaubt, die Basis in Stufe 2 zu verlassen. Die Regel 

nach der getauscht wird, ist dieselbe wie in Stufe 1 . Das 

Optimum ist erreicht, wenn al le relativen Kostenfaktoren 
nichtpositiv C<=0) sind. 

Es laesst sich beweisen, dass in der optimalen Loesung m 
von den x— Spalten sein duerfen. Dazu kann man das Tableau aus 
Fig. 5.1 heranziehen und betrachtet es als pr i ma I es Problem. 
Multipliziert man jede Zeile i Ci = 1 , . . . , rO des Pr i ma I en in 
Fig. 5.1 mit der korrespondierenden Un bekann ten M. des Dualen, 
summiert diese Ze i I en auf und subtrahiert sie \^an der 
Z i e I f unk t i on , so erhae I t man bei einer optimalen Loesung 



1 1 

(o - jZh a n ) *; 

i = 1 

n 

+ (0 * 5Z Mn } x i 



i =1 



(-1 ♦ ir, )v; 



1 ' 1 



(-1 - IT, ) v; ♦ 



n 

(0 - Y~ IT. a. )x* 

v L — i im m 

i = l 
n 

(0 + y TT. a. )x" 



i= 1 



. + (-i ♦ ITJ v; 
•+ (-1 - TT ) v; 







C5.25 



Auf der rechten Seite dieser Gleichung steht die Differenz der 
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Z i e I f unk t i onswer t e der pr imaien und dualen Loesung . Ist die 

Loesung opt i ma I , so s i nd - nach Kap i te I 3.3 — bo i de 

Z i e I funk t i onswer te gleich. Daher die Nu I I auf der rechten 
Se i te . 

Aus den Kuhn-Tucker Rege I n und dem Dua I en kann ersehen 
werden, dass folgende Beziehungen gelten muessen . 



- Hir, a,j * 

i= 1 

tlT-, Qy £ o 



i= 1 



Hir, a =o j = 1 , , 



m C5.3) 



-TT,- >- 

IT; *- 



1 ♦ ITj ^0 
i - \ < 



■1 SIT. *\ 



C5. 4> 



Bei der optimalen Loesung muss also (5.2) nu I I ergeben . 
G I i eder 



Di 



( - Lf. a )x* 



i = 1 



J |J J 



bzw 



(o ♦ Eir, Qij )x; 

i = l j 



verschwinden aber auf jeden Fall, da ja C5 . 3) er f ue I I t ist. Es 
koennen daher alle m x * Coder x~^ groesser nul I, damit 

den 
Ueberlegungen aus Kapitel 4.1 nicht > sein) . 



J J 

Basisvar i abl e, sein Cbeide zugleich koennen nach 
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Alle rest I ichen Gl ieder 

(-1 ♦ I, ) v; 

( _ ! _ TT. ) v- 

ergeben Werte <= Cwegen v? , vT > = und Beziehung CS. 4)). 
Um C5.2) zu erhol len, muessen die Gl ieder einzeln verschwinden, 
es muss daher gel ten: 

V* > ( V.>0 ) I: = 1 

v" > ( v.<0 ) H. = -1 

1 I > 

v: v: =0 m,i * i 

Diese Beziehungen zwischen den v. und IT. werden spaeter mit 
Hi I f e anderer Lieber I egungen auch noch gefunden und dann 
erlaeutert werden. 

An Hand einer Einze Ipunkte i nscha Itung soll dieser 
Algorithmus illustriert werden. 

Gegeben sind die Koordinaten von drei Festpunkten Fl, F2 y F3 
und die Naeher ungskoor d i na t en fuer einen E i nscha I t punk t NP ' 



y Lm2 

-1270.00 

-670.00 

-550.00 

-695. 00 



Weiters sind Strecken und orientierte Richtungen vom Neupunkt 
zu den drei Festpunkten gemessen: 





x CmD 


Fl 


2200.00 


F2 


2645.00 


F3 


2250.00 


NP 


2320.00 
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Strecke CmH Richtung CNeugradD 
NP-FI 587.397 286.9015 

NP-F2 325.955 4.8883 

NP-F3 16 1.013 128.6317 

Mit den bekannten I inear isierten Boobach tungsg I e i chungen 
erhae I t man "Folgendes inkonsistentes, ueber bes t i mm t es 
G I e i chungssys tem Cm i t Dimension in cm und Neuminuten) 

1 0.978910 Ay + 0.204294 ax = 0.9 

2 -0.076697 Ay - 0.997054 Ax - -0.5 

3 -0.900552 Ay + 0.434749 A x = 0.1 

4 0.022142 Ay - 0.106096 Ax = -0.05 

5 -0.194731 Ay + 0.014998 Ax - 0.08 

6 0.17 1893 Ay + 0.356065 ax = -0.09 

An diesen inkonsistenten Gleichungen u/erden, wie ueb I ich, 
Verbesserungen angebracht. Diese Gleichungen koennen dann in 
einem dreidimensionalen Koordinatensystem als Ebenen 
dargestel It werden. In Fig. 5.4 sind die Sehn i t t ger aden 
C I , 2, 3, 4, 5, 6) der Ebenen mit der xy— Ebene eingezeichnet. Es 
koennen nun wie in Kapitel 4.2 die Absolutwerte der Ebenen in 
v— Richtung geb i Idet und summiert werden. Man erhaelt dann 
einen konvexen Polyeder, der an einer Ecke die Optimal I oesung 
besitzt. Auch aus der Theor i e der I i near en Op t i m i erung f o I g t s 
dass die Menge al I er zu I aess i gen Loesungen einen konvexen 
Polyeder aufspannt, und dass die Zielfunktion in mindestens 
einer Ecke ihr Optimum annimmt, 



Als erste zu I aess i ge Anfangs I oesung kann v" = |. 
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Cbzw. v* = |l.| wenn I- negativ ist) genommen werden 



Das Anfangs tab I eau lautet: 



S CNP-FO 
S CNP-F2} 
S CNP-F3) 
R CNP-F13 
R CNP-F2) 
R CNP-F3) 



Bas i s 


L 


Ay' 


ix* 


v i 


0.98 


0.978910 


0.204294 


v 2 


0.50 


0.076697 


. 997054*** 


v 3 


0.10 


-0.900552 


0.434749* 


v 4 


0.05 


-0.022142 


0. 106096** 


V S 


0.08 


-0. 194731 


0.014979 


v 6 


0.09 


-0. 171893 


-0.356065 




1 .72 


-0.23371 1 


1 . 401 107 



Dieser zu I aess i gen Anfangs I oesung entspricht in Fig. 5.4 die 
Stel Ig I. 



Bei der ersten Iteration in Stufe 1 wird Ax + in die Basis 

gebracht (wegen des groessten relativen Kostenfaktors 1 .401 107). 

Das normale Simplexpivot (0.434749*0 f uehr t zu einer 

Approximation, in der Ax* « . I 0/0 . 434749, v* - v^ = 0; das 

heisst, die 3-te Gleichung wird konsistent CI / ). Erhoeht man 

AX + um den Betrag 0.10/0.434749 Cd.h. Ax + - 0.20/0.434749), so 

kann der relative Kostenfaktor von Ax* nochmals verringert 

werden, aber dadurch wird v 3 positiv (man springt vom negativen 

auf den gespiegel t en Punkt im positiven Teil der Ebene 3, dor t 

ist ja der v-Abstand positiv). Das heisst, man ersetzt v~ in 

der Basis durch v* . Dasselbe kann durch folgende 

3 

Ze i lenoperation erreicht werden. Man subtrahiert zwei mal die 
Pivotreihe von der Reihe der relativen Kostenfaktoren. Damit 
wird der relative Kostenfaktor von v" nu I I . Dann wird das 
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Vorzeichen der Pivolre ihe geaendert und v" durch v* in der 
Basis ersetzt. Der relative Kostenfaktor von Ax* ist nun 
0.531609, das heisst, man kann Ax + nochmals vergroessern. Die 
momentane Loesung wird durch I 7/ r epr aesen t i er t . 

Das zweite Pivot C0 . 106096**) entspricht einer 
Approximation mit Ax 1 = 0.05/0.106096, v* = 0, v ^ = 0. In Fig. 
5.4 ist diese Approximation durch I' dargestellt. Bei 
Anwendung der vorhin erwaehnten Ze i I enoper a t i onen reduziert 
sich der relative Kostenfaktor von Ax + auf 0.319417. Daher 
kann Ax* nochmals vergroesser t werden. 

Das dritte Pivot C . 997054*** ) entspricht einer 
Approximat i on, bei der Ax + = 0. 50/0. 997054, v * und v~ nu II s i nd . 
Durch die Zei lenoperation w i r d aber jetzt der Kostenfaktor 
von Ax + -1.674691. Es wird nun an diesem Element C . 997054***5 
pivotiert und v* in der Basis durch Ax + Cbzw . Ax + in der 
Nichtbasis durch v~ ) ersetzt. Die P i vo t e I emen t suche in der 
Spalte erfolgt nach den allgemeinen Regeln aus Kapitel 3.2. 
Die Pivot ierung geschieht mit der sogenannten Rech t ecksr ege I 
Cz.B. NEUMANN K. C19753. Bd. 1. Seite 51,525. 

Man beachte, dass durch das Pivotieren die beiden negativ 
gesetzten rechten Seiten L der dritten und vierten Zeile wieder 
positiv werden muessen . Dies kann man sich leicht an Hand der 
P i vo t suchr ege I und der Rechtecksregel ueber I egen . 
Das Tab I eau nach dem Pivotieren hat folgende Gestalt: 
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Bas i s 


L 


Ay* 


+ 

v 2 


v i 


0.797551 


0.963195 


-0.204898 


Ax " 


0.501477 


0.076924 


1 .002955 


s 


0. 11 80 17 


0.933896«« 


0.436034 


v ; 


0.003205 


0.030303« 


0. 106409 


V 5 


0.072488 


-0. 195880 


-0.015023 


v ; 


0.268558 


-0. 1 44584 


0.3571 17 




1 .259820 


1 .587179 


-0.320361 



Die zu I aess j ge Loesung, die dieses Tableau darstel lt, ist i den t 
mit II in Fig. 5.4. 



In der zweiten Iteration der Stufe 1 wird das Pivot in der 
vierten Ze i le der Spalte von Ay* gefunden. Dies entspricht 
einer Approximation, in der Ay* te 0.003205/0.030303 und v* =* 
und v~ = ist und die durch die Ecke II / dargestellt ist. Der 
Z i e I f unk t i onskoef f i z i en t von Ay* ist nun 1 .526573. Daher 
kann Ay * noch vergroessert werden. Man wandert also 
(vergleiche Fig. 5.4) entlang der Geraden 2 und vergroessert 
dabei Ay + . Das naechste Pivot C0. 933896*0 findet man in der 
dritten Zeile. Hier ist Ay + = 0.118017/0.933896. An dieser 
Stel le kann der Z i e I funk t i onskoef f i z i en t von Ay + nicht mehr 
verkleinert werden, ohne dabei negativ zu werden. Das hei ssl y 
man pivotiert an dieser Stelle und bringt Ay + an Stelle von v* 
in die Basis. Das Tableau erhaelt die Gestalt: 
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Bas i s 


L 


V 3 


+ 

V 2 


v i 


0.675831 


-1 .031373 


-0.654612 


AX* 


0.491756 


-0.082369 


0.967039 


Ay + 


0. 126371 


1 .070783 


0.466898 




0.000624 


0.032448 


-0.092261 


V 5 


0.097242 


0.209745 


0.076433 


+ 

V 6 


0.286829 


0. 154818 


0.424623 




t .060495 


-1 .634628 


-1 .244933 



Die relativen Kostenfaktoren von v + und v* sind —1 .634628 bzw. 

3 2 

— \ .245933, das heisst w iederum / die relativen Kostenfaktoren der 
un t er drueck ten Vektoren v" „ v" sind -0.365372 bzw -0.754067. 
Man beachte, dass die Summe der relativen Kostenfaktoren von v + 
und v" ~2 sein muss . 

Da al le N i ch t bas i svek tor en relatiVe Kostenfaktoren <= haben,, 
ist keine Notwendigkeit fuer Stufe 2 y und das obige Tab I eau 
stellt die optimale Loesung dar Centspricht der Ecke III in 
Fig. 5.4) : 

Ax = . 492 cm 

A y = 0.1 26 cm 

v, = -0 . 676 cm 

v 2 - 

v 3 - 



0.00! Neuminuten 



v 5 - -0.097 Neuminuten 
Vg » 0.287 Neuminuten 



Die Absolutsumme der Verbesserungen ist i .060495. Man beachte, 
dass in der optimalen Loesung v und v nu II sind, das he i ss t y 
die zweite und dritte Gleichung werden konsistent gemacht. 



Dv;1 




FIG. 5.4 
MEDIAN ( 2 - DIMENSIONAL ) 
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Dies aeussert sich in der graphischen Darstel I ung der Fig. 5.4 
dadurch, dass sich die optimal© Loesung im Schnittpunkt der 
beiden Geraden 2 und 3 befindet. 

Wie man aus Fig. 5.4 ersehen kann, ist es der Vor te i I 
dieses A I gor i thmus, dass man mit einer Pivotierung mehrere 
Simplexschritte ueberspringt Cbzw. mehrere Ecken des konvexen 
Polyeders). Eine Pivotierung war nur fuer II und III notendig. 



Zusammenfassend kann ueber diesen Algorithmus gesagt 
wer den : 

a. Stufe I besteht aus m Sehr i t ten y und Stufe 2 beinhaltet 
alle f o I genden Schritte bis zur optimalen Loesung . 

b. Waehrend Stufe 1 wird der Vektor fuer die Basis nur 
aus den Vektoren x + und x" gesucht. Der Vektor, der Nichtbasis 
wird, entstammt nur den Vektoren von v* und v~ . Die Summe der 
relativen Kostenfaktoren von jedem Paar x* und x" ist nu I I, die 
der v * und v~ sind —2. Daher sind diese Kosten jederzeit 
verfuegbar. Sollte in Stufe 1 kein Vektor zu finden sein ; der 
die Basis ver laesst, so ist der Rang der Matrix A kleiner als 
m. In diesem Fa I I wird der Vektor x + Coder xO in der 
Nichtbasis gelassen, keine Pivotierung dur chgef uehr t und fuer 
alle weiteren Schritte ignoriert. Linear abhaengige Spalten in 
A bringen daher keinerlei Probleme; ein grosser Vortei I dieses 
Verfahr ens . 

c. Alle weiteren Sehr itte bis zum Optimum sind in 
Stufe 2. Hier ist nur mehr ein Austausch von Basis v + , v~ mit 
Nichtbasis v * , v~ erlaubt. In Stufe 2 ist eine optimale 
Loesung garantiert, das heisst aber auch, ^in Vektor, der die 
Basis verlaesst, muss immer vorhanden sein. Durch 



AI 



Rundungsfehler kann es al I erd i ngs vorkommen, dass kein Vektor 
gefunden wird. Der Algorithmus ist dann abzubrechen, 
allerdings ist die Loesung meist sehr nahe der optimalen 
Loesung . 

d. Eine wichtige Schlussfolgerung aus diesem Algorithmus 
ist jene. Hat die Matrix den Rang m Cd.h. die Anzahl der 
Parameter x},so sind exakt m der Verbesserungen 
Nichtbasisvariable; also null. Die restli chen Verbesserungen 
sind in der Basis und >*= 0. Daraus folgt insbesondere, dass 
so viele Verbesserungen nu I I sind, wie zur Bestimmung der m 
Unbekannten x gebraucht werden. Oder auch anders ausgedrueckl; 
es sind m Gleichungen konsistent, die restl i chen n — rn 
G I e i chungen be i nha I ten W i dersprueche . 



Dieser Schlussfolgerung laesst sich auch eine interessante 
geodae t i sehe Deutung geben. Diese geodaetische Deutung und ein 
Verfahren, dass auch die schwache Besetzung der A Matrix 
ausnuelzt, werden im naechsten Kapitel behandelt. 
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6. ANWENDUNG DER I - NORM BEI AUSGLEICHUNG VON 
GEODAETISCHEN NETZEN 



6. ! Beobachtungen gleicher Genau i gke i t 



Das in Kapitel 5 vorgestellte Verfahren von Barrodale und 
Roberts kann direkt auf die Ausgleichung von geodae t i sehen 
Netzen unter Minimierung der Absolutsumme der Verbesserungen 
angewandt werden. In diesem Fa I I besteht die Matrix A (n * nrO 
aus den Koef f i z i enten der I i near i s I er ten 

Beobachtungsgleichungen von Strecken, Richtungen u.s.w. Die 
Unbekann t en t e i I menge x besteht aus den Koordi natenzusch I aegen 
der Neupunkte und den Zusch I aegen zu den 
Orientierungsunbekannten. Weiters kann man voraussetzen, dass 
bei geodaelischen Netzen keine linear abhaeng igen Spalten 
auf t r e ten . 

In der optimalen Basis (Loesung) werden daher alle m 
Spa I ten der Unbekann t en t e i I menge x und n — m Spa I ten der 
Unbekann t ente i I menge v sein. Genau m der Verbesserungen v 
werden exakt nu I I sein. Auf ein geodae t i sches Netz uebertragen 
heisst dies, dass m der Beobachtungen keine Verbesserung 
erhalten. Es wird das Netz aus genau diesen m Beobachtungen 
aufgebaut werden. Eine optimale Loesung steht also im 
Zusammenhang mit einer Auswahl von Beobachtungen^ die ausreicht, 
das Netz widerspruchsfrei zu bestimmen. 
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Ein Basiswechsel in der ersten Stufe entspricht dem 
Hinzufuegen einer Beobach t urjg > um das Netz vorerst aufzubauen. 
Ein Basiswechsel in der zweiten Stufe entspricht einem 
Austausch zweier Beobachtungen. Beobachtungen werden dann 
solange ausgetauscht, bis die optimale Netzkonfiguration 
gefunden wurde. Wie schon aus Kapitel 3 bekannt ist, erhae I t 
man jede Bas i s I oesung auch durch Auf I oesen eines m mal m 
Gleichungssystem, gebildet aus den Spalten der Basisvari ab I en . 
Angenommen, es ist eine solche Basismatrix gegeben und 
ausserdem soll die Loesung davon optimal sein. Dann hat man 
ein Gleichungssystem der Form, wie in Fig. 6.1 angegeben, 
auf zu I oesen . 
sortiert, dass eine Einheitsmatrix entsteht) 



CD i e v s i nd der Anschau I i chke i t wegen so 

B 






CD 



AA 







Ml 






I 1 


@ 


<m 





BASIS 



o 



NB 



ZIELFUNKTION 



6 . t 



B 



asismatr ix tuer eine Das i s I oesung 



B< 



I 
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Es ist hier keine Vereinbarung wegen der Vorzeichen von x 
notwendig. Die Toi Imatrizen F und G bi I den zusammen die 



lal 



Ma tr i x A 



e x . 



J 



Cj - 1 , . . . y m) sind in der Basis 



Cwegen der Voraussetzung von I inear unabhaeng i gen Spalten) und 
n — m der Verbesserungen v. Die einzige Vorzeichenvereinbarung, 
die zu treffen ist, betrifft diese v Je nachdem , ob v 
positiv bzw. negativ ist, steht in der dazugehoer igen Spalte — 1 
bzw. + t . Die Indizes B und NB beziehen sich auf die Spalten 
der Verbesserungen, die in der Basis CB3 sind, oder Nichtbasis 
CNB) sind. 

Ein Teil der n Multiplikatoren IT (und zwar n-m) kann 
sofort bestimmt werden. Sie sind ja die Loesung der 
G I e i chungen : 



v B <0 \ 



also sind die Multiplikatoren TT R zt 
Basisverbesserungen bekannt. 



j eder 



n — m 



v B >0 
v < 

B 



Hb = 1 



Da eine optimale Loesung angenommen ist, muessen die 

relativen Kostenfaktoren fuer die Basis nu I I und fuer die 

Nichtbasis <= sein. Daher muss fuer eine v* Spalte gelten: 

NB 

C = " 1 + IL," 
NB 



5! 



bzw. fuer die korrespondierende Nichlbasis v Spalte 



c.- 1 -V„s 



Daraus folgt insbösondere, dass fuer diese Nichtbasis 



NB 



1 ~ K«~ 

NB 



sein muss^ um beide c Cfuer v und v ) <= zu lassen 

NB NB 



Die Multiplikatoren II koennen ebenfalls gedeu ted werden 



Es gibt fuer jede Beobachtung 



bzw 



uer 



jede Verbesserung 



v. C i = t , . . . , n) eine Groesse TT. . Fuer eine Verbesserung 

i j 

v. > ist TT. = I , fuer eine Verbesserung v. < ist TT ~ — 1 

i i i » 

und fuer eine Verbesserung v. = (das he isst fuer eine 
Messgr oesse, die das optimale Netz aufbaut) ist | TI" - f < S5S ! . 

Fuer Verbesserungen v. 7* koennen die Multiplikator en IT • auch 

als Einheitsvektoren in Richtung der Verbesserung gedeuted 

werden. Fuer Verbesserungen v . = kann dieser Vektor jede 

1 

bei i eb i ge Richtung annehmen. Sein Absolutwert ist <= 1 . 

Das in Fig. 6.1 angegebene optimale Basisglei chungssy s t em 
kann nun in folgender Weise geloest werden. 

a. Berechne die Unbekannten x aus.* 



G x = I 



NB 



x = G l 



NB 



C6. 1 ) 



b . setze d i 



1 ese x ein in 



Fx - v n = l 



C6.22) 



52 



und bestimme daraus 



v B = Fx-l B C6.3> 



Ein aehnl icher Wog wird zur Berechnung der Multiplikatoren 
beschr i t t en . 

a. Die Te i I menge IT ist gegeben durch die schon bekannten 

D 

Verbesserungen v . 

B 



TL = 1 wenn v > 

( k = 1 , • ■ -, n - m 

TL =-1 wenn v < 



b. die rest I i chen TT <- 1 <= TT <- I ) , j = 1 , 

NB • NB ■ * 

werden erhal ten aus 



^"mB^X ^ " N B=- (G ' "B «•« 



Mit Hilfe dieses Verfahrens kann auch bewiesen werden, 
dass die in Kapitel I erwaehnten Medianeigenschaften hier in 
einer allgemeinen Form erhal ten bleiben. 

Es sei wieder die optimale Basis aus Fig. 6.1 
vorausgesetzt. Um eine Loesung optimal zu halten, duerfen sich 
die Vorzeichen der relativen Kostenfaktoren nicht aendern . 
Diese relativen Kostenfaktoren koennen bekannt I ich berechnet 
werden aus den Gleichungen*. 
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c-*>1 . 




(-)l 



IT, 



NB 



F D 



Fig. 6.2 

Gl ei chungssystem zur Berechnung der relativen 
Kos tenfak Loren 

Damit c gleich bleibt, darf sich Il R ni ch t aendern, damit 

bleibt auch U gleich und folglich der Vektor c. Nun s t e I I t 

sich die Frage, wie darf man den Beobachtungsvektor I aendern, 

um die TT Q gleich zu belassen und damit M MD und c^ also die 

□ NB 

Optimal i taet zu erhalten. 

Die Antwort ersieht man leicht aus dem Gleichungssystem 



der Fig. 



TT 



bleiben gleich., wenn sich das 



Vorzeichen der Verbesserungen v nicht aender t . CDa fuer TL 
nur das Vorzeichen relevant ist). 



Man kann die Beobachtungen aus der Gruppe l R ver aender n 
Cohne das Vorzeichen der v dabei zu aendern!), die 
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Optimalitaet der Loesung bleibt erhalten. Der Vektor x aendert 
sich dadurch ebenfal Is nicht. 

Wird diese Aenderung als Ausreisser oder als grober Fehler 
interpretiert., dann wird er in jedem Fall erkannt . Natuer I ich 
kann man auch al le Beobachtungen der Gruppe l R in diesem Sinne 
veraendern, die optimale Loesung bleibt erhalten. Es ist 
leicht einzusehen, dass sich dabei nur die Absolutsumme der 
Verbesserungen ver gr oesser t . 

Es sind dies dieselben E i genschaf t en, die den Median bei 
einer direkten Beobachtungsreihe auszeichnen und seine 
Robustheit gegenueber Ausre i ssern bewirkt. Diese Merkmale 
treten hier nur in veral Igemeinerter Form auf. Das heissl, die 
Loesung des Vektors x entspricht einem veral Igemeinerten 
Median. Oder man koennte auch sagen, jene Beobachtungen, die 
das optimale Netz aufbauen, sind vergleichbar einem 
veral Igemeinerten Median. 

Jede Variation — e i ner oder mehrerer Coder auch a] I er) 
nicht zur Med i anb i I düng verwendeter Beobachtungen - wenn man 
dabei auch das Vorzeichen der korrespondierenden Verbesserung 
nicht umformt, bewirkt keine Aenderung der Optimalitaet. 



In einem geodae t i sehen Netz wuerde das bedeuten, dass die 
Zusch I aege zu den Koordinaten und den Orientierungsunbekannten 
gleich bleiben, trotz geaender ter Beobachtungen. Die 
Verbesserungen der geaender ten Beobachtungen erhalten dann neue 
Werte abhaeng i g von den Aenderungen der Beobachtungen; alle 
anderen Verbesserungen behalten ihre Werte. 
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Diese Betrachtungen zeigen die ausseror den t I i che 

Robustheit der Ausgleichung nach der 1— Norm gegenueber groben 
Feh I er n auf . 

6.2 Beobachtungen mit ungleichen Gewichten 



Der Fa I I g I e i chgew i chleter Beobachtungen kann leicht auf 
das Problem ungleicher Beobachtungen gef uehr t werden. 
Hu I i i p I i z i er t man jede der Gleichungen 



1 1 1 



mit dem Gewicht p. der Beobachtung I- , erhae 1 t man 



5~a tj P; Xj -P| V; =Ps l 



\~-\ 



Dadurch wird ein Ersatzproblem mit neuen Beobachtungen 



f. = p. I. ( i= 1 / • -,n) 



neuen Koeffizient« 



a.. = p. a.. 



und neuen Verbesserungen 



( i = ],- ■ -,n ) 
( j = 1,- -,m) 



v. =p i v i ( i = !,• • .n) 



gebildet. Die Minimumsforderung lautet dann 



n n 



Pi V i 



V. I = Min 



i =\ 



i * I 
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Es besteht daher die Aufgabe, folgendes Model I zu loesen: 



ÄX - V = 1 

n 

Y~ | V | = Min 



Das analoge Tableau, zu dem in Fig. 5.1 angegebenen, hat dann 
d i e Form 



i 



^ 


=z^\ 


-0 


D 



o 



o 



= Max 



Fig. 6.3 
Tab I eau fuer gewichtete Beobachtungen 



Mit Hilfe analoger Ueberlegungen koennen die Multiplikatoren JTj 
Ci = 1 , . . . , n) bestimmt werden . Fuer Beobach tungen, d i e 
Verbesserungen v. ungleich erhalten Cfuer jene Spalten, die 
in der Basis sind}, ist der korrespondierende Multiplikator 
TT. ( v * ) = p bzw. TT. (V) = — p . Fuer alle and&r Gn Beobachtungen gilt: 

- P, * IT., * P 



»der 



i i " i 



M*p. 
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Es ist I e i ch t einzusehen, dass die Multiplikator en i m 
Prinzip die Rollo der Gewichte uebernehmen. Fuer 
Beobachlungen, die eine Verbesserung ungleich nu I I erhalten, 
ist der dazugehoerige Multiplikator ident mit dem Gewicht. 
Fuer Beobachtungen, die keine Verbesserung bekommen, ist der 
Absolutbetrag des Mu I t i p I ikators kleiner oder gleich dem 
korrespondierenden Gewicht. 
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ANWENDUNG DER LU ZERLEGUNG BEIM ALGORITHMUS 
VON BARRODALE UND ROBERTS 



D i e Auf I oesung e i nes G I e i chungssys t ems mit H i I f e der 
LU Zerlegung laesst sich beim Algorithmus von Barrodale und 
Roberts vorteilhaft anwenden. Schon BARTELS, R.H. und 
G . H . GOLUB C I 968!) haben d i e Dre i eckszer I egung be i der 
Simplexmethode angewandt. Ihre Untersuchungen ergaben 
guenstige Eigenschaften in bezug auf die Rundungsfehler. 

Es ist wiederum das in Fig. 5.1 vorhandene Prob I ein 
vorgegeben. Generei I sind die Dimensionen al I er Vektoren und 
Matrizen gleich angenommen; also die Anzahl der Beobachtungen 
ist n und die Anzahl der Unbekannten x ist m. 

In Stufe I des Verfahrens von Barrodale und Roberts werden 
nur Unbekannte der Menge x in die Basis gebracht. Es wird nun 
der k~ te Schritt der I . Stufe betrachtet. Dann sind genau k 
Unbekannte x und n — k Verbesserungen v* Cv* ) in der Basis. Die 
Matrix, bestehend aus den Spalten der im k~ten Schritt 
vorhandenen Bas i svar i ab I en , und die rechte Seite (a I so der 
Beobachtungsvektor 13 ergeben das in Fig. 7. 1 dargestel I te 
Gleichungssystem. Dieses Gleichungssystem ist nach dem in 
Kapitel 6.> mit Formeln C 6 . 1 3 bis C 6 . 4) angegebenen Verfahren 
zu berechnen. 
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Die in Fig. 7.1 verwendeten Indizes beim I-Vektor B und NB 
beziehen sich wieder auf jene Beobachtungen, die das Netz 
aufbauen CNB — Nichtbasis) und jene Beobachtungen die 
Verbesserungen erhalten (B~Bas i s) . 



n- k 



ZEILE ZU V--5 





(-.1. 

r-)i 


D° 




@ 



t2 r 



'NB 



. 



Fig. 7.1 
Gleichungssystem der Bas i svr i ab I en in Stufe I 

k- ter Sehr i t t 



Im Ck+1}— ten Schritt kommt die Unbekannte x* in die Basis 
und die Verbesserung v* aus der Basis. Durch diesen 
Basiswechsel wird F um eine Ze i le verkleinert und um eine 
Spalte Cdie der Unbekannten x") ver gr oesser t . Die Matrix G 
wird um eine Spalte und eine Zeile vergroessert. I wird um 
ein Element verringert und l.._ um ein Element erweitert. Das 

N D 

Gleichungssystem mit den neuen Basisvariablen des Ck+1)— ten 
Sehr ittes ist in Fig. 7.2 dargestellt. 
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ZEILE ZU V -- 



k+1 



v^ 



vmutittw 



a\ 



(n- k)-1 




n 



c-)\ 



k+\ 



—da 



L N6 



Fig. 7.2 
Gleichungssystem der Basisvariablen im 
Ck-H )-ten Sehr i tt 

Wie man aus Fig. 7.2 ersehen kann, ist im neuen Schritt k+ t 
zur Berechnung der neuen Unbekannten x im wesent I i chen Gin 
Gleichungssystem der Groesse Ck+O mal (k+l ) auf zu I oesen . 

Allgemein kann gesagt wer den, dass im k-ten Sehr i 1 1 der 
Stufe I zur Berechnung der k Unbekannten x immer ein 
Gleichungssystem der Groesse k mal k zu I oesen ist. Bei jedem 
nachfolgendem Schritt erhoeht sich die Dimension der Matrix G 
um eins. Es kommt immer nur eine neue Ze i le und eine neue 
Spalte zur vorhergehenden Matrix G dazu. Hat man im Schritt k 
zur Auf I oesung von 



Gx = 1 NE 



C7. I ) 



G i n der For 



m 



G = LU 




U = 




zer I eg 



legt und damit (7.13 al 



LUx = l 



NB 



dargestel \t y so kann x in zwei Schritten durch 

zuruecksubs t i t u i er en der beiden Dr e i ecksma tr i zen L und U 
berechnet werden. 



Li aus Lz = l NB mit z = Ux folgt z 
2.j aus Ux = z folgt x 



Bei der Dre i eckszer I egung koennte man das in WILKINSON J.H. 
C 1 965) auf Se i te 226 angegebene Verfahren m i t Ze i I enlausch 
verwenden. Es ist numerisch stabi I und fuer Rechenan I agen 
besser geeignet als herkoemml i che Methoden. Die Elemente der 
beiden Dre i ecksma tr i zen koennen am selben Platz gespeichert 
werden wie in der Matrix G Cunter Ver nach I aess i gung der 
D i agona I e I emen t e C—l) von L ) . 
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Es ist also im k-ten Schritt die LU Zerlegung der Matrix G 
vorhanden. Im Ck+O-ten Schritt waechst G um eine Zg i le und 
e i ne Spa I te . Fuer d i e Neub i I düng der LU Zer I egung von G i m 
Schritt Ck+1) brauchen dann nur mehr die Elemente der neu 
hinzugekommenen Spalte und Ze i le berechnet zu werden. 




12ZZZZZZZZZZZZZZZZZZZZZZ 



'IßltmiUUULL 



k + 1 



NEU ZU BERECHNEN 



Fig. 7.3 
Neu zu berechnende Elemente von G = LU 
bei Sehr i tt Ck+1 D 

Waehrend Stufe 1 des Bar r oda I e ' sehen Algorithmus waechst 
die Matrix G von der Groesse I mal 3 auf m mal m an . Nach 
genau m Schritten ist das Ende der Stufe 1 erreicht und Stufe 2 
beg i nn t . 



In der Stufe 2 werden nur mehr Basis— und 
N i ch t bas i sverbesserungen ausgetauscht. Die Matrix G hat 
waehrend der ganzen Stufe 2 die gleichbleibende Dimension 
m ma I m . 
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(t)1 




V 







!A 



NB 



-k lj 



Fig. 7.4 
Situation in Stufe 2, alle m der Unbekannten x 
und n-m der Verbesserun g q n in der Basis 



Kommt waehrond eines Schrittes in Stufe 2 die Verb 



esserung v 



in die Basis und v A aus der Basis, so bewirkt dies einen 
Ze i lenauslausch im Gleichungssystem der Fig. 7.4. Wieder kann 
man \for\ der Voraussetzung ausgehen, dass die LU Zerlegung von G 
des vorhergehenden Schrittes verfuegbar ist. Tauscht man in 
der Matrix G* eine Zeile aus, so braucht auch hier nicht die 
ganze Dr e i eckszer I egung gemacht zu werden. Es genuegt, die in 
Fig. 7.5 angegebenen Te i le der schon LU zerlegten Matrix neu zu 
berechnen . 
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AUSGETAUSCHTE 
ZEILE 



NEU ZU 
BERECHNEN 



Fig. 7.5 
Neu zu berechnende Te i le der LU Zer legung 
von G waehrend Stufe 2 

Natuerl ich kann die jewei is vorhandene Dr e i eckszer I egung 
von G auch dazu benuetzt Uerden, die Multiplikator en zu 

ber echnen . 



G T TT 



NB 



F T TL 



C7.23 



Mi t 



wird C7 .2} zi 



LU 



(LU)' TT, 



NB 



U T L T 



NB 



F'rt E 

F T IT. 
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und IT kann in zwei Schritten berechnet werden. 
NB 

1.) aus U T z = - F T IT B mit z = L T IT NB folgt z 
2) aus L T TT NB = z folgt T NB 

Fuer die Suche nach der Verbesserung, die aus der Basis f ae I I t, 
wird in beiden Stufen des Algorithmus von Barrodale und Roberts 
das Produkt der inversen Basismatrix mit der urspruengl i eben 
Spalte (jener, die in die Basis kommt!) benoetigt. Dieses 
Produkt entspricht dem Vektor der neuen Bas i s var i ab I en in den 
Tab I eaus des Simplexalgorithmus und auch im Tab I eau beim 
Barr oda I e ' sehen Algorithmus. 

Bekannt ist die Basismatrix, geb i Idet aus den Spalten der 
dazugehoor i gen Bas i svar i ab I en . 



-)\ 







+C— } I sol I wieder andeuten, dass bei Verbesserungen v > hier 
- I und bei Verbesserungen v < hier +1 steht. 
Man hat aus 



:*-)! F 



G 





/ 1 
ä B 


_ 


a B 1 




1 ä NB 

i 




a J 



C7.35 



den gesuchten Vektor 
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a 



Qnb 



zu berechnen. 

Durch Ausmultiplizieren von C7.3) wird erhalten 



i!iIa B + F a NB = a B 



Ga NB = a NB 



Die dreieckszerlegte Matrix von G ist immer verfuegbar, es kann 
daher aus 



LUä =p 

NB *NB 



der Te i Ivektor a berechnet werden. Durch Einsetzen in 

NB 



i!i I a B = a 8 - F a NB 



wird a und somit der gesamte gesuchte Vektor erhalten. 



Es koennen noch verschiedene Vereinfachungen auftreten. 
Der urspruengl i che Vektor, bestehend aus den Te i Ivektoren a Q und 

a.._,ist eine Spalte aus der Matrix der Beobach tungsg I e i chungen / 

NB 

solange man sich in Stufe 1 befindet. In dieser Spalte treten 

viele Nu I lelemente auf. Wenn zum Beispiel a aus lauter 

NB 

Nu I lelementen besteht, dann gi I t 
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(*-)I a B = a B 



In Stufe 2 werden nur mehr Spalten der Verbesserungen zum 

Bas i swechse I verwendet . In diesem Falle ist a _ ein Vek t or m i t 

b 

Nullelementen. a ist ein Vektor mit Nu I lelementen. und an der 

NB 

korrespondierenden Stel le der gewechselten Verbesserung, die in 
die Basis kommt^ steht eine 1 . 
Dann gilt 



LU ä 



NB 



Ml ä B = -Fa N8 



Die Suche nach den Variablen., die dem Basiswechsel unter I iegen, 
geht ana I og nach den Rege In, wie sie in der 
Standardsimplexmethode angewandt werden. 

Die Matrix G ist bei geodae t i sehen Netzen in der Regel 
schwach besetzt. Der Vor te i I der Dr e i eckszer I egung besteht 
darin,, dass auch die LU zerlegte Matrix von G weiterhin schwach 
besetzt bleibt Cim Gegensatz zur Inversen von G ,, die in der 
Regel dicht oder fast vo I Istaendig besetzt sein wird). 



Die Matrix G wird am besten komprimiert abgespeichert 
C Vandom packing") . Dabei werden nur die Ni ch tnu I I e I emen t e der 
Matrix in einem Vektor Cin bei i eb i ger Reihenfolge) und in zwei 
Hi Ifsvektor en die jeweiligen Zeilen- und Spaltenindizes des 
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Elementes abgespeichert 



•I 3 



O 



o 



SPALTENINDEX |S *©*©2O©&0## 



ELEMENTE G • • • • 12>3 • •©••} (#000000 

ZEILENtNDEX |Z 



> 


• 





















• 000000 



Fig. 7.6 

Vandom pack i ng" einer schwach besetzten Matrix 

Wie in Fig. 7.6 angedeu ted, steht also das Element 12.3 der 
Matrix G in Zeile 3 und Spalte 20. Der Rest der drei Vektoren 
G y IZ y IS ist mit Nullen aufgefuellt. 

Bei dieser Art von Absp'e i eher ung lassen sich leicht 
Elemente h i nzuf uegen . Sie brauchen nur am Ende der 
Nichtnul lelemente von G y IS,IZ angehaengt zu werden. Bei der LU 
Zerlegung muessen sehr o f t auch Ze i I en getauscht werden. Dazu 
brauchen nur die Indizes an den entsprechenden Stel I en in IZ 
und IS gewechselt zu werden. So I I ein Element als geloescht 
gekennzeichnet werden, so ist es guenst ig, die entsprechenden 
Elemente in den Vektoren IZ und IS negativ zu setzen. An 
diesen Plaetzen koennen dann wieder neue Elemente eingetragen 
werden . 



Die LU Zerlegung besteht hauptsaechl ich aus inneren 

Produkten ^or\ Zeilen— und Spaltenvektoren. Es muss also in 

derse I ben We i se auf Spa I ten und Ze i I en zurueckgegr i f fen werden 

koennen. Dies ist bei diesem Abspe i eher ungspr i nz i p 
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platzsparender moeg I ich^ im Gegensatz zu Abspe i eher ungsar t en 
Cz.B. mit verbundenen Listen), die einen Ze i I en— oder 
Spa I t en i nd i zesvek tor und einen Poi ntervektor Cder zum naechsten 
Element der Zeile oder Spalte weist) haben . Bei der Bildung 
der inneren Produkte werden dann saemtliche Mu Itiplikationen 
mit nu I I unber uecks i ch t i g t bleiben. 



Durch diese Art von Abspeicherung wird also die schwache 
Besetzung der Matrix G und die bei der LU Zerlegung bleibende 
schwache Besetzung ausgenuetzt. Ausserdem werden saemtl i che 
Operationen mit Nullelementen unUrdrueckl . 
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8. PUNKTTRANSFORMATION UNTER MINIMIERUNG DER 
ABSOLUTSUMME DER RESTKLAFFUNGEN 



8.1 Problemstellung und mathematische Formuli er ung 

Gegeben ist ein System von n Punkten P. und ein System von 

n Punkten 0. Auf das System der Punkte Q. so I I eine 

Transformation ausgeuebt werden, deren Parameter so zu 

bestimmen s i nd, dass die Absolutsumme der resultierenden 

Abstaende (Restkl affungen) zwischen den P. und Q. ein Minimum 

wird. Die Transformation der Punkte Q. so I I im folgenden aus 

i 

einer Mas s t absaender ung , Verdrehung und Verschiebung in der 
Ebene bestehen. Man beschraenkt sich also bei der Darstel I ung 
auf die sogenannte He I mer t tr dnsf orma t i on in der Ebene. 
Aenderungen und Veral I gerne i nerungen des Transformationsgesetzes 
bereiten jedoch keine Schwierigkeiten. 



Die Koordinaten \/on P seien x ,y , jene von Q seien vor 

iii i 

der Transformation gegeben durch x. +dx. y y. +dy. und nach der 
Transformation durch x. +v . , y. +v . Weit er s wird angenommen, 

dass d i e Groessen dx , , dy. , v . , v . klein gegenueber den 

Koordinatendifferenzen der Punkte P. seien. Das heisst, die 

Punkte 0. sind schon durch eine vor I aeuf i ge Transformation mit 
den P. zur Deckung gebracht worden. 
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Bezeichnet d^> den Ver dr ehungsw i nke I , dyu den Masstabsfaktor 
und de j de d i e Versen i ebungen i n x— bzw . y — R i chtung, so gilt: 

x y 



x;+ v xi = ( Xj+ dxj) cosdf - ( yj + dy ; ) sin dip + x ; du + dc ) 
y; + v yi = (Xj + dX;) sin dip + ( y- t +dyj) cosdy + y t du + de, 



und in erster Naeher 



ung 



v,i = x.d^ -y-.dip + dc x + dx, 
v y i = yjdu +Xjdip + dc y +d> 



i= 1 ,--,n 



C8. IcO 



Als Res tk I af fung wird definiert 



p, =A < * \ 2 



<i yi 



p >0 



C8. Ib) 



und es wird gefordert 1 



n 

) p ; = Min 



ui 



C8. Ic) 



Diese Optimierungsaufgabe C8.1a) > C8 . 1 b ) , C8 . I c) kann, wie man 
unmittelbar einsieht, durch folgende gleichwertige ersetzt 
wer den 

v xi = Xjd^A - yjdip +dc x + dXj 
v yi = y|d^u + Xjdcp +dc y +dy ; 



C8.2a) 



^ v 2 + v 2 ' < 

V I VI 



xi y i n 



j> =M 



in 



1 = 1 



C8.2b) 



C8.2c) 
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Gleichungen (8.2a), (8.2b) ergeben zusammen mit der 
gewuenschten Minimumsbedingung (8.2c) eine n i ch t I inoarö 
Optimierungsaufgabe. Die Relation in (8.2b) besagt nur, dass 
man den transformiertön Punkt innerhalb oder an der Peripher ig 
eines Kreises mit dem Radius p zulaesst. 

Diese so f ormu I ierte Optimierungsaufgabe (8.2a), (8.2b), 
(8.2c) kann aber wiederum durch eine I ineare 
Optimierungsaufgabe ersetzt werden. AI lerdings hat dann 
letztere unend I ich viele Nebenbedingungen. Eine Folgerung 
davon ist, dass die urspruengl i che Optimierungsaufgabe C 8 . 2a ) , 
(8. 2b), (8.2c) jedenfalls konvex ist. Der Ueber gang von 
(8.2a), (8. 2b), (8.2c) zu einer li near en Op t i m i erungsauf gäbe 
mit unendlich vielen Nebenbedingungen geschieht 
f o I gendermassen . Fuer v . und v 1 . gilt 



» x i ~ PiCOS Ci Richtungswinkel der 



V, = p-, sin^ 



t; 



Restklaffung 



durch Multiplizieren mit cosÄ bzw. sinA, wobei A den ganzen 
Bereich von bis 2JT ueber s tr e i ch t , 



V xj COSX = pjCOST, COSÄ 

* A < 2 TT 



v yi sinX = PiSinfj sin Ä 
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jnd durch Addieren wird erha I ten 



v xi cosX + v yi sin Ä = p-, cos(^ - Ä ) 



bzw 



v xi cosX + v yi sin X < pj £ X < 2 F 



C8.3) 



Das heisst, man ersetzt 



x i vi — r i 



durch 



v xj cos A + v yi sinA ^ Pi ^ A < 2 TT 

Deutet man v . , v. vor ueber gehend als Koordinaten eines 
Punktes, so bedeutet C 8 . 1 b 2) y dass der Punkt auf der Peripherie 
eines Kreises mit dem Radius p. liegen muss. Glei chung C 8 . 2b) 
hingegen verlangt nur, dass der Punkt im Inneren oder auf dem 
Umfang I i egen muss, Gleichung (8.2b) ist demnach aequivalent 
dem System unend I ich vieler Ungleichungen (8.3), also 



v xj cosX + v ;SinX < p ^ X < 2 IT 



C8. 4) 



Jede dieser Ungleichungen (8.4) spannt eine Halbebene auf, die 
von einer Tangente an den Kreis begrenzt wird und den Ursprung 
enthaelt. Der Winkel £" ist die Nor ma I r i ch tung der Tangente an 
den Kr eis. 
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Fig. 8. I 

Darste I I ung der Halbebene fuer die Nor ma I r i ch tung tf 

Die Zwangsbedingungen C8. 1b) werden durch unend I ich viele 
Ungleichungen C8.4) ersetzt. Das heisst geometrisch^ dass der 
Kr e i s durch den Durchschn i t t von unend I i ch viel en 
Tangen t enha I bebenen ersetzt wird. Setzt man fuer v . und v 



y i 



gemaoss (8 . 1 a) i n den Ung I e i chungen C8 . 4) e i n y so erhao I t man 



(X|C0SÄ +yj sin X) du + (XjSinX ~y i cosA)dip+ cosX dc x * sinA dc y + 
+ dxi cos X + dyi sin X ^ p; 

bzw. = picosf^i -X) = pi für ri = A 
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oder auch 

( X; cos X + y ; sin X ) d u, + ( x, sin X - y, cos X ) dip + cos X dc x + sin X de , 
- Pi ^ - ( dXj cos X + d yj sin X ) 

( i= 1, •■ • ,n) 

< X < 21 



C8.6) 



m 



i t 



Y P» - Min 

i= 1 
Werden nun fuer jeden Punkt nur endli ch viele X .. gewaehlt 
Cj = \j . . .,m. 5 j so ergibt sich 



(xiCOsXy + y i sinX ij )du + ( x ; sin X^ -y.cos \ V} ) dip + dc x cos X l; +dc y sin \ V} 

-pi- -(dXjCOsXij +dyiSinXij) 



1 1 
^p, = Min Pi > 



i =i 



C8.7) 

1 ■ • ,n j= 1 ,-■ -, rrii 



In C8.7) wird der Kreis durch ein Tangentenpolygon 

approximiert, und die X.. sind die Normo I r i ch t ungen an die 

»J. 

Pol ygonse i ten. 

Aufgabe C8.7) stellt aber auch ein li near es Programm dar, 
welches das urspruengl iche nichtl ineare Programm (8.1), (8.2), 
approximiert. Wird Aufgabe C8.7D als PRIMALE I ineare 
Optimierungsaufgabe betrachtet, so kann auch das dazugehoer i ge 
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DUALE Problem erhalten werden. 
Es I autöt : 



n m; 



- ZIIZ ( x i cosX ü + y-. sin h ) TTij 

n mj 

- IZ ZI ( x i cos x ü - yi sin a ü > h 

■. -. i j = i 

n m*, 

- H5Z TijCosXy 

n i"i 



H i 



m - 



i = i 



i) 



».. 






1=1 



n m; 



= 

= 





< 1 

< 1 



H 51. ^ dX ' C0S X Ü + d Vi S ' n X Ü ) ^Ü = M 



ax 



F.. > 



C8.85 
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8.2 Anwendung des Satzes vom komp I emen t aer en Schlupf 

Der Satz vom komp I emen t aeren Schlupf besagt, dass fuer 
jene Beziehungen des Pr i ma I en (8.7) mit der Relation < die 
korrespondierende Unbekannte des dualen Problems nu I I ist. In 
der folgenden Diskussion wird angenommen,, dass in der optimalen 
Loesung bei p. > die zu den optimalen Richtungen £*• 



9Q 



hoerenden A-. C = t- ) — Spa I ten vertreten sind. Dies wird bei der 



spaeter zu beschreibenden Iterationsverfahren mit hinreichender 
Genau i gke i t erzwungen . 

Befindet sich eine Variable des primalen Problems in der 
Basis Cund zwar handelt es sich um die Basis einer optimalen 
Loesung), so ist der Schlupf der entsprechenden Ungleichung des 
dualen Problems gleich nu I I (und umgekehrt, da ja das Duale 
wieder als Primales und das Primale wieder als Duales 
aufgefasst werden kann) 



Es ergeben sich folgende Moeg I i chke i ten: 
a . ) Punkt mit P ( > Cp ( . . . Bas i svar i ab I e). 

Ist ein p | > 0, so kann fuer diesen Punkt i nur eine der 
Gleichungen (8.7) mit = erfuellt sein Cnaemlich fuer A • ■ » — £*; ) . 
Alle anderen G I e i chungen muessen mit < erf ue I I t se in. AI so 
tr i t t in d i esen G I e i chungen Seh I upf auf . Demnach s i nd alle IT.. 
nu I I C f uer alle j ungleich j*). Da die Summe der Tl.. den Wert 
I ergeben muss, folgt notwendigerweise TT-, ~ 1 . 

Dies gilt eigentlich streng nur fuer das Problem (8.6). 
Bei dem endlichen Problem (8.7) kann die Relation = auch fuer 
zwe i Bed i ngungen e i ntreten Cd i es i st dann der Fa I I , wenn der 
Punkt an eine Ecke des Tangentenpolygons zu I i egen kommt). In 
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dem spaeter zu beschreibenden Iterationsverfahren wird solange 
i teriert, bis sich die beiden Tangentenrichtungen genuegend 
genau naehern . Geometrisch bedeuted d i es, dass der Kreis mit 
Radius P| in der Richtung X-, durch zwei (annaehernd) gleiche 
Tangenten ersetzt wird. Die Summe der zwei auftretenden TT. 



muss 



se in. Sie k 



oennen 



dah< 



er als ei nz i ges 



TL 



■J 



t d 



er 



Groesse 1 aufgefasst werden. Der Grund dafuer ist, dass man 
ein nichtl ineares Problem durch ein I ineares ersetzt hat. 



b. ) Punkt mit pj = 0. 

Die Unbekannte p des Pr ima len ist in der optimalen Basis, 
er hae I t aber den Wer t nu II. Die dazugehoer i ge Seh I upfvar i ab I e 
des Dualen muss daher verschwinden. Es koennen mehrere TT* 
auftreten. Ihre Summe muss eins sein. 

Die Unbekannte p. des Pr ima len ist nicht in der optimalen 

Basis Csie ist Ni chtbas i svar iab I e), hat also ebenf a I Is den Wert 

nu II. Die korrespond i erende Seh I upfvar i ab I e kann > se i n„ 

muss aber nicht. Daraus folgt, dass die Summe der hier 

auftretenden TT < 1 sein kann. 

»J 

Die zu den auftretenden (1 gehoerenden R i chtunqen X. 
sind verschieden. In einigen Fae I len kann es vorkommen, dass 



h i er nur ein Tl.. der G 



r oesse 



1 Cp. ...Basis) auftritt oder d 



ass 



die Schlupfvariable ist Cp- ...Nichtbasis). Dies tritt bei 

Transformationen auf, wo mehr als zwei Punkte die Restklaffung 

bekommen. In der Regel werden aber bei Pj ™ 

CO. . . . N i chtbas i s !) zwei verschiedene Richtungen X.. , die Summe 

der Tl.. kleiner 1 und die Seh I upfvar i ab I e > sein. Bei 0=0 
ij ^ r i 

CO.. . .Basis) treten dann drei verschiedene Richtungen X- auf 

r, «» ij 

und die Summe der dazugehoer i gen IT.. ist 1 . Geometrisch 

»J 

bedeuted dies,, dass dann der Punkt durch drei Tangenten 
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e i ngezwaengt ist. 

Zusammenfassend gilt: 

Punkt mit 0. > TI .= 1 , IT.. « . . . j * J • 

» 1 ij <J 

k 

Punkt mit 0; = CO- . . .Bas i svar i ab I e) IZTT • - 1 

» 1 r • . — - »j 

Punkt mit D- = ( 0. . . . N i chtbas i svar i ab I e) 5~ TT-- <= I 
r, r, ^ u 

geht ueber jene k TT.. > 0, d 



CD 



er oummat 1 ons 1 ndex 



1 e i n e 1 nem 



Punkt mit p = auftreten) 



Es koennen die Unbekannten IT.. des Dualen als Laengen von 

Vektoren mi t Richtung A-. aufgefasst werden. Fuer einen Punkt 

mit D. > ist der Vektor p Cm i t . den Komponenten TT.cosX..., 
r i > ij* «j 

TT- sin X..- bzw. Ip.l = O in Fig. 8.2 dargestellt. 




gcnte 



Fig 8.2 

Unbekannte des Dualen gedeuted als Vektor 
fuer einen Punkt P i > 



In einem Punkt mit p = koennen die auftretenden IT.. auch zu 
einem einzigen Vektor zusammengef ass t werden. Fig. 8.3 zeigt 
die Situation fuer einen solchen Punkt mit drei Tangenten. 
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3 Tangenten 




Fig 8.3 
Punkt mit der Restklaffung null 
drei Tl.. zu einem Vektor zusammengef ass t 



Das sichere Kriteri um fuer das Auftreten der Restklaffung nul I 
ist dann gegeben, wenn |p|< 1 . 



8.3 Deutung des dualen Problems 



Werden wie in Kapitel 8.2 alle Unbekann ten TT- des Dualen 
als Vektoren p. aufgefasst, dann gilt 



fuer pi > = 



I 



Pi 



TT: COSf; 



F. sin V. 



PJ = i 



f uer p ; = : 



' u;, cos t-, I 



T n Sinti, 



i 



I i2 cosr l2 
TT i2 sinr, 2 



H i3 cosf i3 
Ti 3 sinf i3 



I] COS fj 

TT, sin fj 



|Pd=s 1 



und das duale Prob lern C8.7) kann in folgender Weise 
interpretiert werden: 

Suche Vektoren p. y so dass gilt: 



I . 5 



-E 



i--i 






p-, = o 



die Summe der Inprodukte ist nu I I 



2.) 



-L 



r=l i J \ 



P; = 



die Summe a I I er Exprodukte ist nu I I 



C8.9) 



C8. io: 



3.) 



-51 Pi =0 

i = 1 



C8. i n 



die Summe der Vektoren p. ist nu I I 
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4. ) 



|P;I± 1 



C8. 12) 



Was bedeuten nun diese vier S 



ummen 



Man kann sich vorstel len, dass in jedem Punkt Q| eine 
Kraft p. der Groesse |p.| angreift. Es ist dann moeg I i ch, den 
Summen (8.9), (8.10), (8.11), (8.12) eine mechanische Deutung 
des dualen Problems (8.8) zu geben. 



id 3.3 



T P; =o 



i = i* 



Die Summe a I I er angrei fenden Kraefte (resul t ier endo Kraft) muss 
null sein. Das heisst y die Kraefte duerfen dem als starr 
aufgefassten System keine tr ans I a t or i sehe Beschleunigung 
er teil en . 




Fig. 8.4 
angreifende Kraft p 
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id 2. ) 



-E 



Pi =o 



I =1 



Das Exprodukt zweier Vektoren kennzeichnet das Drehmoment. Die 
in Q- angreifenden Kraefte ueben ein Drehmoment aus. Die Summe 
dieser Drehmomente ueber alle Punkte i Ci = I , . . . 4 n ; 
resultierendes Moment) muss nu I I sein. Das heisst, die Kraefte 
duerfen dem System keine Rotationsbeschleunigung ertei I en . 



>d I . ) 



L 

i =i 






.p, = o = 



rix, 



im 



i = 1 




•OP; 



Hier muss angenommen werden, dass das System der Punkte nicht 
mehr ganz starr ist, sondern vielmehr einer Mass t absaender ung 
(Streckung^ Schrumpfung!) keinerlei Widerstand entgegensetzt. 
Die angreifenden Kraefte muessen dann so sein, dass keinerlei 
den Masstab veraendernde Beschleunigung auftritt. 

Eine geodaetische Deutung von C8.9) — C8.12) ergibt sich 
aus der inneren Feh I er t heor i e eines Punkthaufens CVergleiche 



dazu MEISSL P. C 1969)5. E; 



i nd dann d 



I 



Punk t ver seh i ebungen aufzufassen, von denen zu fordern i st, dass 
ihr gemeinsamer Ver seh i ebungs- , Verdrehungs— und Mass t absan t e i I 
g I e i ch null ist. 



Zusammenfassend kann gesagt werden, dass die Punkte Q; 
durch Kraefte p. nach den P- geschoben werden, ohne dass ein 
Ver seh i ebungsan t e i I, ein Ver dr ehungsan t e i I und Mass t absan t e i I 
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uebr j g bleibt - unter der Forderung: minimiere die Summe der 
dabo i erhal tonen Res tk I af f ungen . 



8.4 Loesung und Beschreibung des Iterationsverfahrens 

Das approximierende I ineare Programm C8 . 7 ) loest man am 
besten ueber sein Duales C8.8D mit Hilfe der r ev i d i er t en 
S i mp I exmethode . Es f ae I I t dabe i auch d i e Opt i ma I I oesung von 
C 8 . 7 ) an . Dabei werden fuer jeden Punkt P Ci = 1 , . . . , n } 
einige X- Cj = 1 y . . . , m. ) gewaeh I t . Guens t ig ist es, zu 
Beginn den Kreis z.B. durch vier Tangenten zu approximieren, 
also etwa die vier Richtungen mit ., 100, 200, 300 Neugrad zu 
waeh I en . Man formul i er t ein erstes approximierendes I ineares 
Programm C8.8) und loest dies mit Hilfe der revidierten 
Simplexmethode. Es wird eine erste Loesung fuer die Iteration 

, , ,_(n ,_m ,„(1) ,«(i) = (n sin , , . 

k = 1 du , d^ . dc x , dc y , P^ - • • - - Pn erh a I ten . 



Aus den Transformationsparametern koennen gemae: 



C8. IbD 



(l) 



.n> 



entsprechende v. , v. , bzw. K I af f ungen 



P'^A (v<Y + (v y <;>r 



berechnet werden. 
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Dieses erste approximierende I ineare Programm wird die 
eigentl i che Aufgabe C8.7) noch nicht genau genug angenaeher t 
haben^ daher gilt: 

CVergleiche dazu die geometrisch© Darstell ung in Fig. 8.55. 

Fuer jeden Punkt konstruiert man eine neue Ungleichung 
C 8 . 7 } , bzw. im dualen Problem wird ei ne neue Spa I t e 
h i nzugef uegt . Der zu dieser Spalte gehoerende Richtungswinkel 
A.. # (zu jeder Spalte des Dualen gehoer t ja ein bestimmtes A.. ) 
kann berechnet werden aus 



cos Xjj. = 



sin X ;j . 



ii ) 



YÜ 



,(U 



,i\) 






C8. 135 



Die zu jedem Punkt gehoerenden Indizes m. Csie zeigen ja an, 
wie viele Spalten im ersten Iterati onsschr i t t momen tan fuer 
jeden Punkt vorhanden sind) werden bei jeder Iteration um eins 
erhoeht. Es gi It also 



(k*t) (k) 1 

m = m + l 

i i 



Damit hat man fuer jeden Punkt i eine neue Spalte Cmit X - « ) 



bekommen. Die neuen Ungleichungen werden 
augenb I i ck I i chen Loesung n i cht erf ue I I t . 



von 



der 



86 




Fig. 8. 5 
Situation im Punkt i bei Itorat i onsschr i t t k 
Tangentenpolygon bestehend aus 5 Tangenten 



Geometrisch bedeuted dies Cvergleiche auch F i g . 8 . 5 *J : 

(k ) 
Der transformierte Punkt 0. I iegt in dem zu den Richtungen X - 



Cj 

Inkr e i sr a 



(k 



I , . . . , m . } gehoerigen Tangentenpolygon mit d< 



_(k 



d i us p. > aber nicht im Inkreis selbst. Das Polygon 

wird durch eine weitere Seite mit Normo lr i chtung X... ergaenzt, 

( k ) 
sodass jetzt Q. ausserhalb des neuen Polygons I iegt. Die 



Anzahl der A.. fuer Punkt i ist jetzt m 



ikt!) (k) 

= m . 
i i 



.m 



Dualen werden dazu einige Spalten h i nzugef ueg t . Bei der 

rev i d i er ten S i mp I exmethode kann d i e augenb I i ck I i che Loesung 

._(*) ,_(k) ._(W) ,_(k) -(kl -(kl . . - .. 

da , dtp , de , de , p . . . , p als Ausgangs I oesung t uer d i e 

Berechnung des erweiterten Problems dienen. Damit ist der 

Sehr i t t k abgesch I ossen . 
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In di 



osom o»nno i tarier 



ter 



t man so lange, bis p 



^(k) 



p. un t er eine 



9 



ewisse Schranke f ae I lt„ oder besser: bis sich die aus v, 



(k) 



x \ 



v berechnete Richtunq C der Restklaffung, der zu der 

y i ö i 

Basisspalte ij des Dualen gehoerenden Richtung A. genuegend 
genau naehert. 



8 . 5 Rev i d i er ter S i mp I exa I gor i thmus be I Loesung des Dua I en 



Das Duale Problem hat "Folgende Form (fuer jeden Punkt 
gleich viele A . . , j = I , . . . , m ) . 



Unbekannte des Oualcn IT- 



'J 



•a 



für Ä | 



für h. 



z^\ } 


^\ 2 1 



für \ . 







1 -.. 

j "'□, 

' 1 


1. 

"i 



JL L 



J 



n Schtupfvanable 



-11 



'In 



•21 



'2n 



1 [ 



Tnl 



Max 



Fig. 8.6 
Form des dua I en Prob I ems 
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Eine Spalte (zu jeder Spalte gehoert eine Rieht ung X .. ) 
besteht aus folgenden Elementen: 

- (Xj cos X ;j + y t sin \ Vj ) 
-(Xj sin k H - y. cos X,. ) 

- cos X^ 

- sin X i; 







C8. 143 







dx. cos X.. + dy. sin X.. 



Die Multiplikatoren lauten dann 



m 



T = ( dyu , dip , dc x , dc y , p, , • • , p n ) 



Diese Multiplikatoren werden so bestimmt,, dass gilt 



rwA 



C B = C B - J_ m k Q kB =0 B Basis 



n*4 



} m a Q — NB • Nichtbas 



C NB C NB L — '"h ^kNB 
k = 1 



Die rel at i ven Kostenfaktoren c werden geb i I de t aus 

|J 

Ci = I , . . . , n ; j = 1 „ . . . , m 5 : 
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Cjj = ( dx ; cos Xj. + dy ; sin \ ) - 

- C - ( X; cos X ;j + y ; sin Xj. ) dyu - ( x, sin X^ - cos X t/ ) dip 
- dc x cos Xjj - dc y sin X. (j + p. ] 



Durch die spezi e I ie Anordnung der Spa I ten kann hier eine 
Vereinfachung bei der revidierten Simplexmethode erreicht 
werden. Jede Basismatrix B (sie wird geb i Idet aus den Spalten 
der Basisvariablen) hat bei der Loesung des dualen Problems 
eine spez i e I Ie Form. 



G° 


® 


p i 


G2) 



Fig. 8.7 
Basismatrix B bei der Loesung des Dualen 



F hat die Dimension C4 mal n3, und jede Spalte besteht aus den 
ersten 4 Elementen von C8.14D. F ist demnach ei ne Auswahl \zon 
Spalten der A • Matrizen Cj = 1 , . . . , m ) von Fig. 8.6, wobei die 
Spalten so sortiert sind, dass die 1 Elemente von C8.14) 
zusammengereiht eine Einheitsmatrix ergeben. Aus den 
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Ueber I egungen von Kapitel 8.2 ergibt sich fuer C eine 4 mal 4 
Matrix. Auch die Spalten von C bestehen aus den ersten 4 
E I ementen von C 8 . 1 4 } . Die Matrix D ist schwach besetzt und 
besteht aus einer Auswahl von Spa I ten von I. 

Es kann nun zwischen drei verschiedenen Problemen 
unterschieden werden. 

Prob I em I : 

Zu berechnen sind bei vorgegebenen Basisvariablen Cdas heisst, 

bei bekannter Basismatrix B) die Unbekannten Tl.. des Dualen. 



Es ist ein G I e i chungssys tem der Form 



F C 



I D 



V 



IT 



c 



gegeben. Die Indizes bei II und IT zeigen die Groesse der 
beiden Vektoren an. Durch Ausmultiplizier en er hae 1 t man 



ir„ ♦ 



d ir t = d 

F n d - D TT, 



f ir n + c u, = c 
f ( d - d ir 4 ) + c ir n = c 

Fd + (C- FD) TT 4 = c 
TT = (C-FD)"' (c - Fd) 



und dam i t durch d i e spez i e I I e Form von c und d 
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c = 



/ 1 



(C- FD) F 



F n = d - D T t 



Die Matrix CC — FD) ist eine 4 mal 4 Matrix, deren Inverse 
leicht zu bestimmen ist. 

Prob I öm 2 : 

Zur Bes t i mmung der Spa I t e > die aus der Basis faellt^ benoe t i g t 
man die Vektoren q und t, wie sie in Kapitel 3.4 (revidierte 
Simplexmethode) aufscheinen. q ist bekannt I i ch m i t dem Vektor 
der derzeitigen Bas i s I oesung ident. 



IT = 



IT. 



\ 



= q 



un 



d t wird er ha I t 



en aus 



/f 



c 



D 



/ . I 



1W 



la.l 



a, 



n 

\ l 
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dadurch ergibt sich 



t t ~- (C - FD)" (a, - Fa n ) 



t„ = a n "Dt, 



Dabei ist wieder zu beachten, dass a n ein Einheitsvektor ist 
und dadurch noch Vereinfachungen entstehen. 

Die Suche nach der Spa I Le, die in die Basis kommt, erfolgt 
nach den Regeln aus Kapitel 3.4. Ebenso kann zur Berechnung 
der relativen Kostenfaktoren Csie bestimmen ja die Spalte, die 
in die Basis kommt) wieder die strukturierte Form der 

Basismatrix B ausgenutzt werden. Die Unbekannten des Pr i ma I en 
Cd la, d^> , de, de, p , . . . , p } , also die Parameter der 
Transformation, koennen wieder durch Ausnutzen der speziellen 
Form der Basismatrix B vereinfacht berechnet werden. 



'f t 



C T D T 



I \ 



\ p i 



1 \ 


C J 


\ C W 



Hier beinhaltet x die Transformationsparameter und im Vektor p 
sind die Rostk I af fungen p , . . . , p . Die Vektoren c und c. 
beinhalten die Faktoren der Zielfunktion jener Spalten vom 
Dualen, die in der Basis sind. 
Durch Ausmu I t i p I i z i eren w i rd erha I ten 
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J r-T 



= (C T - D'F') (c 4 - D'cJ 
= c - F T x 



CC — D F ) ist aber die Transponierte der vorhin verwendeten 
4 mal 4 Te i I i nver sen CC — FD) . Um d i e Spa Ite zu bestimmen, 
die in die Basis kommt, benoetigt man die relativen 
Kostenfaktoren. Sie werden berechnet aus 



C = 



F T I 

C T D T 



\ / \ 
x 



a l so aus 



c n = c n - F' x - p 
c 4 = c 4 - C T x - D T p 



Wie aus dem bisher Gesagten leicht zu sehen ist, benoetigt 
man eigentl ich immer nur die Inverse von CC — FDD bzw. deren 
Transponierte. Diese Matrix uebernimmt praktisch die Ro I le der 
Inversen von B. Damit stellt sich die Forderung nach einer 



vereinfachten Berechnung von 
Bas i swechse I . 



CC - FD) 



nach e i ner 



Prob I em 3 : 
Se i jetzt 
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l a « 



die neue Spalte, die \n die Basis kommt 
Elementen von C8.14). 



besteht aus den 



Fal ! a: 
Eine der Spalten von C, D unterliegt dem Basi saust ausch . Dabei 
bleiben F und I unveraender t . Die neuen Tei Imatrizen koennen 
durch eine Modifikation dargestel It werden. Seien C, D die vor 
dem Basiswechsel und C, D die nach dem Basiswechsel vorhandenen 
Matrizen und j der Spa I t en tausch i ndex innerhalb von C Calso 
entweder 1, 2, 3 oder 4), so gilt: 



C = C ♦ (a 4 - C.j } ej 



e . C4 ma I 

i 



I D Einheitsvektor mit 1 an der 
Stel le j des Basiswechsels 
j — te Spa I te von C 



D kann leicht modifiziert werden, es braucht nur die I in der 
entsprechenden Spa I te getauscht zu werden. Dies ist aber auch 
nicht notwendig, wenn die beiden getauschten Spalten demselben 
Punkt angehoeren. Dasselbe kann auch durch folgende Operation 
erreicht werden. 



D = D +(e k -e, ) e? 
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j ist wieder der Spaltenindex innerhalb der Matrix C y al 



so 



entweder I , 2 , 3 oder 4. e (in diesem Fall ident mit a ) und 
e. sind zwei Einheitsvektoren der Groesse n mit 1 an der Stelle 
k bzw. \ y wobei die 1 an der Stel le I dem Punkt entspricht, dem 
die neue N i ch t bas i sspa I t e angehoer t . Entstammen beide Spalten 
demselben Punkt, so sind k und I ident und D bleibt 
unver aender t . Gesucht ist dann 



[C ♦ (a, -C, ) ej - F(D + (e. -e, ) e. T ] 

** J J Kl 



-1 



bzw . ausmu I t i p I i z i er t 



[C - FD +1 (a 4 -C. -F. k ♦ F., )e]l 



-1 



M ( 4-4 ) 



x ( 4 • I 



y T ( 1 * 4 ) 



Dazu verwendet man die vereinfachte Formel fuer 



~] 



(M + xy T ) = M 



-i 



1 +x r M _1 y 



M xy M 



C8. 155 



D i e Ma t r i x 



M"' -- (C - FD ) 



ist dabei schon aus dem vorherigen Simplexschritt bekannt 



Fal I b: 
Eine der Spa Iten von F, I wird ausgetauscht . In jedem Fall 
wird dabei F geaendert . So I 1 ten die auszutauschenden Spa I ten 
demselben Punkt angehoeren, so bleibt I unveraendert und auch C 
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un 



d D. Die Matrix F wird modifiziert in der Fori 



F + ( a 4 - F.j ) e T j 



Hier ist wieder F.. die j~te Spalte Calso jene, die aus der 

Basis f ae I I t ) und e. ein Einheitsvektor mit an der Stelle j . 

J J 

Gesucht wird wieder 



fC - (F - (a t - F., ) e ■ ) D ]"' 



bzw . 



[ C - FD ♦ (a 4 - F.j ) e. D ] 

[ i i . i , 1 



M ( U * ii ) 



14*1) y'(1-4 ) 



und es kann wieder die neue Inverse von CC — FD) mi t Forme 
C8. 15) berechnet werden. 



Gehoeren die beiden Spalten., die dem Basiswechsel 
unter I i egen, n i cht dem g I e i chen Transformat i onspunkt an, so 
aender t sich F, und die Matrix I muss umsortiert werden, um 
w i eder E i nhe i tsmatr i x zu se i n Cdazu muss aber e i ne Spa I te aus 
C, D mit einer von F., I getauscht werden). Dann hat man wieder 
die Form der Basismatrix von Fig. 8.7 erreicht. Durch diesen 
Tausch aender t sich aber auch noch C, D, und es muessen F, C und 
D in der bisher besprochenen Art modifiziert werden. In diesem 
Fal I wird sich keine grosse Vereinfachung ergeben und es wird 
besser sein, die Inverse von CC - FD) gleich direkt zu 
berechnen . 
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Bemerkung : 

Im Prinzip uebernimmt die Inverse von CC — FD) die Ro I le 
der i nversen Bas i smatr ix B . Da man alle Spa I ten des Prob I ems 
immer neu generieren kann,, ergibt sich ein Speicherplatzbedarf, 
der mit der Punktanzahl waechst. 



8.6 Geometrische Deutung der relativen Kostenfaktoren 



Die relativen Kostenfaktoren c. haben die Fori 

»J 



i dx. cos X.. + dy. sin X.. + (x. cos X.. + y. sin X.. du + i 

r i ' » 

I + (x s sin X;. - y ; cos X 4j ) d<p" + de* cos X Vl + dej sin \ V] \ - p* = c 



Lz pj cos( V\ - X,j ) 



dyu*^ d^p*, c ' c x* " dc y * , p|* sind die Loesung einer 

S i mp I ex i t er a t i on y die noch nicht optimal i s t , waehrend pj die 

aus diesen Tr ansf or ma t i onsgr oessen errechnete Restklaffung ist 

Cmit Gleichung C8.1) berechnete v. , v und mit Restklaffung 



P. = \K' 



werden . 



+■ v 



y 



). Die Situation kann geometrisch dargestel It 
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Fa I I a: Nicht optimale Loesung 




A:q einer Spalte , die sich 
in der Basis befindet 



T< 



c. >n ( nicht in der Basis ] 
}NB 



Fig. 8.8 
Fa I I einer noch nicht optimalen Loesung 



Wenn die Loesung noch nicht optimal i st, liegt der Punk t Q . * in 

Fig. 8.8 noch ausserhalb des Kreises mit dem Mittelpunkt in Pj 

und dem Radius p * . Die Richtung ?? . der Restklaffung faellt 

nicht zusammen mit der Richtung A. D , die zu der Spalte gehoer t , 

i b 

die sich in der Basis befindet. Im S i mp I ex t ab I eau bedeutet 

dies, dass c. > ist. 
iNB 
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Fal I b: Situation bei gefundener Optimal I oesung 



Xjg( Basis ] 



/ 



\ 




X- ( Nichtbasris) 

c- < 
l NB 



Fig. 8.9 
Opt i rna I I oesung 

Bei der optimalen Loesung fallen die Richtung der Restklaffung 
Cv • 4 v • } und d i e zur Spa I te der Bas i svar i ab I en IT . gehoerende 
Richtung A. R zusammen. 



8.7 Aendern der Ausgangsk I af fungen dx, dy eines Punktes 



Eine optimale Loesung aender t sich nicht Csiehe auch 
Kapitel 3), wenn saemtl i che c.. ihre Werte beibehalten. 



?enauer y es muss t uer 



di e Bc 



— und fuer die Nichtbasis 



c .. <= ö erhalten bleiben 

»J 



00 



Aender t man die Ausgangsdaten dx, dy in einem Punkt auf 
dx + Adx y dy -+* Ady, so bekommen die relativen Kostenfaktoren 
eine neue Form. 

c f uer e i ne Bas i svar i ab I e : 

der ursprueng liehe relative Kostenfaktor ist 

c = ( dx cos V + dy sinr ) + 

+ ( x cos 'fr +y sinf ) du + ( x sin?^ - cosf ) dy + dc x cos£- + dc y sin 2: 
- p- 

f uehr t man anstatt dx — dx + Adx, bzw. dy— dy + A dy ein, so 
ergibt sich ein neues c* 

c* = (dx cos tr + dy sin ^ ) + ( Adx ebs^ + a dy sin?" ) + 

+ (x cos^+ y sin^ ) dju Wx sin^ -ycos'M dw + de cos£- + de sin^ 
-P - Ap = c 

damit nun 5 = 5* bleibt, muss sein 

Adx cos tr + a dy sin tr - a p =0 

d i es i st dann der Fa I I , wenn 



Adx = Ap COS f 

Ady = Ap sin v 



Analog dazu darf sich das Vorzeichen von c Cj . . . f uer al le 



101 



N i chtbas i s i nd i zes) im Punkt nicht aend 



ern 



c- - (dx cos X + dy sin \- ) + 

+ (x cos Xj + y sin )\- ) du + (x sin l- - y cos Aj )dip fdqcosAj +dc y sin A 
- p < 



das neue c" m i t den Groessen dx + Adx, dy + Ady erg i bt 



c- = ( dx cos Äj + dy sin Aj ) + ( a dx cos Xj + a dy sin Äj ) + 

+ (x cos Aj + y sin ^Jdyu+fx sin A r y cos Aj ) dip *dc x cos X ] + de, sin Aj 
- p - a p < 

Damit bewiesen werden kann., dass auch hier c* <= bleibt, 
ersetzt man den Ausdruck 



Adx cos Xj * a dy sin Xj - Ap 



i ndem man 

AdX = Ap COS V 

Ady = a p sin ^ 

einsetzt und erhae I t dann 

Ap cos£r cos Xj + Ap sin^ sin Xj - a p = a p oos(f - X. ) 

damit ergibt sich fuer c * folgende Form 



c = p cos ( v - X. ) + A p cos ( f - X. ) - p - a p ^ 

^(P + Apjcosff-Xj) -p - Ap <0 

= (p ♦ A P )(cos(f -Xj) - 1 ) ^ 
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das heisst, c" wird ebenf a I Is <- bleiben. 

Zusammenfassend kann gefolgert werden: 
Wenn man einen Punkt Q- (mit ß. > 0) in seiner f Richtung vom 
Punkt P. entfernt, so aendern sich die Unbekannten du, dip, dc x , 
dc y der Transformation nicht. Die Restklaffung in diesem Punkt 
steigt um den Betrag der Verschiebung 

Dies ist das Ana I ogon zu den Eigenschaften des Med i ans bei 
direkten Beobachtungen und bei dem im Kapitel 6 behandelten 
Fa I I des Ansatzes der verm i ttelnden Beobachtungen. 



Es kann der Fa I I auftreten, dass die Loesung des Dualen 
nicht eindeutig wird. Und zwar wenn ein p- =0. Fuer diesen 
Punkt s i nd saemt I i che c = Cf uer Bas i s und auch fuer d i e 
Nichtbasisvariablen). Wuerde man eine N i ch t bas i svar i ab I e in 
die Basis bringen., so aenderte man dadurch den Wert der 
Ziel funkt i on n i cht . Man kann be I i eb i ge II. fuer d i esen Punkt 
in die Basis bringen, die Optimali t ae t geh t dadurch nicht 
verloren. Auch die Multiplikatoren Cdas heisst, die 
Unbekannten des Pr imalen) bleiben konstant. Wie man leicht 
sieht, muss dies der Fall sein, denn fuer einen Punkt mit der 
Restklaffung p. = muessen die Gleichungen (8.7) des Pr i ma I en 



uer a 



lle V. die Relation = haben 



Grundsaetzl ich kann gesagt werden, dass bei einem 
Auftreten eines Punktes mit p. = ein veral I gerne inerter Median 
gefunden wurde (also etwa vergleichbar mit dem Medianwert bei 
einer hessreihe mit ungerader Anzahl von Beobach t ungen ) . 
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8.8 Bedingungen fuer Loesbarke i t des dualen Problems 



Interpret i er t 



man saem 



t I i che Loesungen II.. > des Dual en 



als Vektoren p. , wobei gi It 



o 



TT cos f. 
I sin e*. 



a. 



Pil = l = 1 



jnd 



0: 

/ 



Pi = 



TT., cos t;, 



T^sinc",, 



1 

U. 2 cos ?. 2 






+ 


H i2 sinr i2 j 





V 05 ^ 



Tf i3 sinc- 13 



I cosc~ 



IT., sinr 



a. 



i i 



*\Ä2Äl <1 



Pil = TTi *1 



dann koennen die ersten vier Beziehungen des dualen Problems 
wie folgt geschrieben werden. 



ET I, cos t\ = 

i = 1 

t TTi sin^i = 

i = t 

E TTj X; cos^ + TT-, y.sin^ 
£ -ITj y, cosf, + TT i x i sin£ 



1=1 








Loest man dieses Gleichungssystem nach zwei bei i eb i gen Punkten 
01 und 02 auf, erhae I t man 
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Hn cos V \z cost c 



02 



IT, sin £; + TT sin t n 

01 Ol 02 02 



x TT cosr + x TT cos?" +y J[ sinZ" * y TT sinf 

01 01 01 02 02 02 7 01 Ol 01 7 02 02 02 



- y TT cos2- - y TT cos^ + x TT sin2r * xtf sin? 



01 01 



0! y 2 02 2 Ol Ol 



01 02 02 



02 



n-2 



cos Z: 



= -51 TT. sint: 

i 

= -H x.TT. cos^. + y TT. sin^ 
= -5~-y.ir.cos£ + x.TT sin2T 

* J i l i i i i 



D i e Summen auf der rechten Se i te gehen ueber d i e rest I i chen n — 2 
Punkte. In Matrizen zerlegt, ergibt dies 




X 01 X 02 Xdi ^02 



'^01 "^02 X 01 X 2 




n-2 



- LI IT-, cos^ 

i = 1 
n-2 

- Y_ \ sin 2Tj 

i = 1 
n-2 

- / x. TT. cos^. +y. TT. sin f. 

i = 1 
n-2 

- 5~" -y I cos?r + x T sinf 

L J i i i i i i 



i = 1 



Legt man nun ein neues Koordinatensystem x, y zugrunde, wobei 
die x— Achse in der Verbindungsl inie von 01 und 02 I iegt und der 
Ursprung sich in der Mitte von der Verbindung 01 und 02 
bef indet, weiters der neue Masstab so gewaeh I t wird, doss die 
Strecke 01 -Ursprung und 02-Ursprung gleich einer Einheit des 
neuen Systems ist, so wird obige Matrix orthogonal, wenn man 
sie zusaetzlich noch normiert. Fuer 01 und 02 gilt dann, bzw . 
fuer das neuQ Koordinatensystem x , y ■ 
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X , = - 1 



X 02 = 



y 01 = ° 
y 02 = ° 



2 - 2 - 2 - 2 



R =V xqi + v 2 + y ; + y 02 



•^ 



Normierungsfaktor 



Das Gleichungssystem erhae I t dann folgende Fori 




"VJ ° 







f2 



VI 



o 



o 




i 

V2 



n-2 

2_ TT; sin^i 



i = 1 

n-2 



VI 

1 
- p ) x.I. cos£-.- y. TT. sin£- 

* ■*- i = 1 



ausmu I t i p I i z i er t 




2 5Z (-1 + x, ) TTi cosZ'i - yiTIi sin?Ti 



n- 2 



Y -( 1 + X;) IT, cos f., + y i lT i sin^i 



i = 1 
n-2 



^ (-1 + Xj ) F. sinr + y- TT- cosr 



i = 1 






i = i 
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und fuer beide Punkte 01 und 02 getrennt 



/ | 



y 01 



cosCj -sin^ 



sinf, cosK ] 



■Ux, 






IT, 



cosr n _ 2 -sinr n . 2 



sin^ n . 2 cosf n . 2 



+ x 



n-2 



y n - 



i. 



n-2 



C8. 165 



02 



>, 



■ oosf sint 



■smC, -cosfr,. 



+ x 



IT 



* I 



cosf , sinf , 

n-2 n-2 



sinr n _ 2 -costr n _ 2 



Wx 



n-2 



Yn-2 



TT 



n-2 



Stel I t man die Frag© nach den Bedingungen, dass in den 
zwei Punkten 0i und 02 p = wird, so muss gel ten 



I D I < 1 I D i < 1 

I K oi I J ro2 1 

Unter der Annahme, dass sich alle rest I i chen in der Basis 
befinden, sind nach den Betrachtungen von Kapitel 8.2 al le 
rest I i chen TT™ I G I e i chungen C8 . 1 6} koennen dann ver e i nf acht 
geschrieben werden. 



2 P 01 = R , U , + R 2 U 2 + • 

2p 02 = r ; u ; * R 2 u ; + • 



+ R u 

n-2 n-2 

■+ R" ,u' , 

n-2 n-2 



C8. 17) 



Der Vektor 2p ergibt sich aus der Summe von um t- gedrehten 

Vektoren u. . Die Vektoren u. sind aber die Verbindungsvektoren 

» i 

vom Punkt 02 zu den resti i chen n-2 Punkten des Systems. 
Analog ergibt sich der Vektor 2p Q aus der Summe von 



um V. + 200 gedrehten Vektoren u! 



u 



ind aber d 



Vektoren vom Punkt 01 zu den restl i chen n-2 Punkten, 



107 



Die nachfolgende Fig. 8. 10 veranschaul icht die Lage der 

Vektoren u. , u ' , der Koordinatensysteme x, y 

Curspruengl iches) und x, y (neues), wei t er s der Drehwinkel 

(Restkl affungsr i chtungen) V- . 




Fig. 8.10 
Lage des neuen Koordinatensystems, der Vektoren u 
und u .' und der Richtungen Ü. . 



!08 
8.9 Bedingungen fuer drei zu transformierende Punkte 

Fuer drei zu transformierende Punkte vereinfachen sich 
Gleichungen (8. 17) zu: 



2P , = R , U , 0,: Pa, = 
2P„ 2 = R;u,' 02: p o2 =0 



P,: P| >0 

Um in den Punkten 01 und 02 p = zu erreichen,, muessen die 

Abso I u t be t r aege I p Q j < 1 und |p 02 | ^ 1 sein. Das wiederum 

bedeuted, dass die gedrehten Vektoren u. und u' kleiner 2 sine 

Dies ist aber nur moeg I ich y wenn der dritte Punkt in einem vc~ 

t 
zwe i Kre i sen Cm i t Rad i us 2 und mit Mittel punkt in 01 bzw . 02^ 

begrenzten Bereich liegt. Fig. 8.11 veranschaulicht die 

S i tua t i on . 



Verallgemeinert heisst es, dass immer die zwei am we i tesi 

auseinander! i egenden Punkte die Restk I af fung p = bekommer 

Also bestimmen diese beiden Punkte die Tr ansf or ma t i onsgr oesse- 
du, dvj? , dc x , dc y . 



Befindet sich ein grob falscher Punkt im Inneren diese- 

beiden Kreise,, so kann er aufgedeckt werden. Liegt er 

ausserhalb, so erhae I t or bei der Transformation p = 0. Se ' 

Fehler wird auf den ihm am naechsten I i egenden Punk 1 
ueber tragen . 
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Fig. 8.11 
f eh 1 er auf deckender Bereich fuer Transformation 
mit drei Punkten 



Bemerkung : 

D i es er Fall ist etwa zu verglei chen mit einer Str eckenmessung . 
Wird die Strecke nur einmal gemessen, so kann keine Aussage 
ueber Fehler gemacht werden. Wird sie zweimal gemessen,, so ist 
es auch nicht moeglich zu unterscheiden, welche davon richtig 
ist. Eine Aussage ueber etwaige grobe Fehler ist erst dann 
moeg I ich, wenn mindestens eine zweifache Ueber bes t i mmung 
vorhanden ist Cdie Strecke also dreimal gemessen wurde). Im 
Fa I I einer ebenen Punk t t r ansf or ma t i on muessen somit mindestens 
vier Punkte vorhanden sein. 
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8.10 Konf i gura tion von Punkten, in der ein Punk t die 
Res tk ! af f ung null erhalten muss 



Nimmt man zum Beispiel eine Anordnung von vier Punkten, in 
der ein Punkt sehr weit entfernt von den anderen ist, dann 
bekommen Gleichungen C8.163 die Form: 



2 P 01 = R , u ^ * R 2 U J 2 

2p 02 = r>; tt, + r;u^it 2 



C8. 183 



Trifft man vorerst noch keine Vereinbarung darueber, in welchem 
Punkt p = se i n so I I , so gilt ob i ge G I e i chung C8 . 1 8) noch ganz 
allgemein. Die einzige Bedingung, die bes teh t y betrifft die 
Absolutbetraege der Vektoren p. Sie muessen in allen vier 
Punkten P. kleiner oder gleich 1 sein. 
A 1 so 

|Pb,l * IP02! > TF , ^2 * ] 

|P.I = IT. JP2M2 



Man nimmt zwei der drei nahe beieinander I legenden Punkte als 
jene Punkte an, die das neue Koordinatensystem x v y bestimmen. 
Es stellt sich nun die Frage, wie weit muss der vierte Punkt 
von den anderen entfernt sein, um auf jeden Fal I die 
Restklaffung p = zu bekommen. Fig. 8.12 soll die Situation 
veranschaul ichen. 



1 1 1 



P 2 ( £ 2 -- ) 



V r 50 9 




Fig. 8. 12 



er zu transformierende Punkte 



Es so I I zuerst nur der Vektor p m fuer den Punkt 01 betrachtet 
werden. Die Gleichung dazu lautet: 

2p 01 =R,u 1 TT, ♦ R 2 u 2 lT 2 (|p ,| < 1 ) 



Wie man aus der Skizze ersieht, ist es gar nicht moegl ich, dass 
in dieser Konfi gur a tion die Bedingung I P 01 1 < = 1 erfuellt wird. 
Nach dem dualen Problem muss aber der zu jedem Punkt i 
Ci = i ,...,, rO gehoerende Vektor |p| < = 1 sein. Sonst waere 
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das duale Problem in der Form C8.8D nicht I oesbar . Egal wie 

die Wahl der Reslk I af f ungsr i ch lungen t- x ist „ es wird nie 

moegl ich sein, 2|p 01 | <= 2 zu erreichen. Die einzige 

Moeglichkeit besteht dar in, IL so klein zu waehlen ( IL < I ) , 

dass 2 1p..) < ä 2 wird. J\ 2 < 1 heisst aber |p j < I und das 
wiederum bedeutet, dass p 2 = wird. 

In dieser Art von Konfiguration wird der Punkt P~ in jedem 

Fa II p 2 = haben . G I e i ches gilt auch f uer den Vektor p 

Vera I I gerne i nerung : 



Es kann auch eine al I gerne ine, von den Richtungen f. der 
Restk I af f ungen unabhaeng i ge Regel fuer das sichere Auftreten 
eines Punktes mit p = gegeben werden. Wie leicht aus den 
vorherigen Abschnitten zu ersehen ist, haengt es von den 
Entfernungen der Punkte untereinander ab. 

Es kann f o I gende a I I gerne i ne Rege I abge I e i tet werden : 
Ist die Differenz der Entfernung von 01 nach einem Punkt P* 
Cvon 02 nach P*) minus der Summe der Entfernungen von 1 nach 
al I en anderen Punkten, ausser nach P* und 02 Cvon 02 nach al I en 
anderen Punkten, ausser . P* und 01 }, groesser als der Abstand 
zwischen 0! und 02, so ist P* ein Punkt, der die Restklaffung 
nu I I erhae I t . 



Diese Regel ist ganz unabhaeng i g von den Richtungen der 
Res tk I äff ungen & 1 '■ und zeigt den Extremfal I auf. Ein Punkt P*, 
der eine bestimmte Entfernung von einer Gruppe anderer Punkte 
hat, kann demnach die Restklaffung nu I I erhalten. Aus der 
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obigen Regel laesst sich leicht ablesen, dass dies umso 
wahrscheinlicher zutrifft, je weiter der Punkt P * von der 
Gruppe der anderen Punkte entfernt ist Cman denke auch an die 
mechanische Deutung des Dualen im Kapitel 8.35, 

Durch eine guenstige Ver te i I ung der 
Res t k I af f ungsr i ch t ungen kann natuerl ich diese extreme Grenze 
noch herabgesetzt werden, denn wie schon gesagt, die erwaehnte 
Regel ber uecks i ch t i g t nicht die C . und bi Idet den Extremfal I . 
We i ters muss noch gesagt werden , dass es schon genuegt, wenn 
diese Regel fuer einen der Punkte 0t oder 02 erfuellt ist. 



8.11 Punkttransformation mit ungleicher Gewichtung 

Natuerl ich kann auch eine Punk t tr ansf or ma t i on unter 
Einfuehrung von Gewichten g. Ci = 1 , . . . , n; n Anzahl der Punkte!) 
dur chgef uehr t werden. 

Das pr i ma I e Problem C8.7) bleibt in seiner Form erhalten 
bis auf die Z i e I f unk t i on . Die Z i e I f unk t i on lautet: 



C 9i Pi = Min 



i= \ 



mit p. >« 

Es kann w i eder das dazugehoer i ge dua I e Prob I em erste I I t werden 



1 I 4 



n m; 



) ) ( x. cos Ä .. 4- y. sin Ä .. ) F. 

ist i«' 



n m 



.= 1 j = 
n m 



i =1 j = 
m 



YZ Z_ (x, sin X.. - y. cosXy ) TT, 

ist j=l 

n nnj 

C LZ TT ti cosA-. 



" IZ JZ "Tu sin X 



r i, 



= 

= 

= 

= 



i=i 



Z nj 



=S 9n 



j"1 



Z Z (dx.cosÄ. + d y. sin X .. ) TT,. = Max 
• = i j= i 

Wendet man hier wie in Kapitel 8.2 den Satz vom komplementaeren 

Seh I upf an, so f o I gen aehn I i che Ergebn i sse und zwar : 



Punkt mit Pj > 






Punk t mit Pj = Cpj ... Basisvar iable) 



Eir,- 



9; 



j=1 



Punkt mit p.= Cpj . . . N i chtbas isvar i ab I e) / TT-. < = 



9 ; 



i = 1 



CDer Summa t i ons i ndex j geht ueber jene k TT.. > 0, die in ein< 
Punkt mit' p. = auftreten). 



Deutet man wieder die TT-, des Dualen als Vektoren p. , so gilt: 



Punkt mit p. = 



Punkt mi t p. > 



IPil * *'. 



!Pil = 9; 



Diese Vektoren p Cbzw . die Unbekannten des Dualen!) 
uobernehmon praktisch die Ro I le der Gewichte in dem Problem. 
Auch die veral I gerne i ner ten Medianeigenschaften bleiben in der 
selben Weise erhalten,, wie sie in Kapitel 8.7 behandelt wurden. 
Eine Verschiebung der Ausgangsk I af f ungen eines Punktes Coder 
mehr er er oder aller Punkte) mit p > entlang seiner t- Richtung 
um einen Wert A p aendert nichts an der Optimal itaet der 
Loesung - die Transformationselemente bleiben gleich und die 
Absolutsumme der Res tk I äff ungen steigt um den verschobenen 
Be t r ag . 

Auch fuer die Loesbarkeit des Dualen koennen aehn I i che 
Bedingungen erhalten werden. Beziehungen C8 . 1 6) wer den dann 

zu ' 



Poi = y ( R, U, IT, + 



R n-2 U n-2 ^-2 ) 



P02 = -r ( r; u; \ + 



+ R' , u' , IT , ) 

n-2 n-2 n-2 



Es ist klar zu sehen, dass man durch geeignete Wahl der 
Gewichte in 01 und 02 immer p = erreichen kann. Es gelten 
w i eder ana I oge Bez i ehungen Ca I I e M. C i = 1 , . . . , n~2) s i nd g I e i ch 
den g. , wenn in 0t und 02 die Restklaffung nu I I sein so I I). 



Po, = y ( R, u , g, + 



R n-2 U n-2 9 n . 2 



P02 = Y { R ' U ' 1 9 ' + 



+ R:.2 U n-2 9n-2 ) 



Um in den Punkten 01 und 02 p =* zu bekommen, 



muss a I so 



Pen < 9 



01 



Po, < 9 



02 



sein. Dies kann entweder durch Wahl von grossen Gewichten in 
01 und 02 oder durch kleine Gewichtsannahmen in den restlichen 
n — 2 Punkten erreicht werden. 
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9. PROGRAMM ZUR AUSGLEICHUNG EINES LAGENETZES NACH DER I -NORM 



Zur Aufstellung des Problems (4.2} aus Kapitel 4 verwendet 
man die bekannten I i near i s i er t en Beobachtungsgleichungen fuer 
Richtungs— und S t r eckenbeobach t ungen, wie sie al I gerne i n beim 
vermittelnden Ausgleich vor liegen. Die Koeffizienten koennen 
analog nach dem vermittelnden Ausgleichsschema bestimmt werden. 
Die rechte Seite ist fuer den linearisierten Fall, wie ueb I i ch , 
gemessene minus gerechnete Beobachtung. Es ist nun moeglich, 
in aehnl icher Weise wie beim vermittelnden Ausgleich, die 
Or ientierungsunbekannten zu "eliminieren' 1 . 



9.1 "Elimination" der Orient ierungsunbekannten 

Man haelt die Koordinaten der Neupunkte einstweilen fest 
Cd.h. die Koordinateninkremente sind null) und orientiert die 
R i ch tungsbuende I Cunter der Forderung- minimiere die 
Absolutsumme der Verbesserungen). Es sind dann fuer jedes 
Buende I d i rekte Beobachtungen vorgegeben und d i e Loesung f uehr t 
auf den Medianwert. Dies wird nun fuer jedes R i ch t ungsbuende I 
dur chgef uehr t . Auch bei der Ausgleichung nach der 2~ Norm gibt 
es eine solche Methode^ nur werden dort die arithmetischen 
Mittelwerte gebildet. 

Im Algorithmus von Barrodale und Roberts laesst sich dies 
folgendermassen erreichen, Man beschraenkt die Wahl der 
Spalten, die in Stufe 1 in die Basis kommen, nur auf die 
Spalten der Orientierungsunbekannten. Bei Anwendung des in 
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Kapitel 6 angegebenen Verfahrens Cmit Formeln (6.1) bis (6.4)) 
ergibt sich die Vereinfachung,, dass die dort verwendete Matrix 
G immer Einheitsmatrix bleibt. 

Es werden also die Medianwerte fuer jedes R i ch tun gsbu ende I 
bestimmt. Damit hat man saemtliche Orientierungsunbekannte in 
der Basis und die dazugehoer igen Verbesserungen In der 

Nichtbasis. Berechnungen Cwie z.B. Bestimmung der Unbekannten, 
Multiplikatoren, relativen Kostenfaktoren usw.) brauchen noch 
nicht gemacht zu werden. Man hat somit etliche Schritte Cso 
viele, wie Orientierungsunbekannte vorhanden sind) fuer Stufe 1 
schon erledigt, ohne dabei umfangreiche Berechnungen 
durchfuehren zu muessen . Erst dann steigt man in den normalen 
Algorithmus ein und bringt die r es t I i chen Unbekann t en 
(Koordi nateninkremente) in die ' Basis, beendet Stufe I und 
beginnt mit Stufe 2. 

Man folgt mit dieser Methode zwar nicht ganz genau dem 
Algorithmus von Barrodale und Roberts, wie numerische 
Untersuchungen spaeter aber gezeigt haben, bleibt die Anzahl 
der Schritte in Stufe 2 in etwa gleich gross. Immerhin bringt 
dies eine Einsparung an Rechenzeit von ca. 30—50%. 



9.2 Programme GENBGL und MAVENA (ASCII FORTRAN) 

a.) GENBGL (Generiere die Beobachtungsgleichungen) 
In diesem Programm werden die Elemente der rechten Sei te, 
der Koeffizientenmatrix und die Naeher ungswer t e der 
Beobachtungen berechnet. Diese Daten werden auf einem 
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D i sc— F i le in komprimierter Reihenfolge abgespeichert. 

b.) MAVENA CMinimiere die Absolutsumme der Verbesserungen 
eines vermi ttelnden Netzausgleiches!) 

Dieses Programm leistet die 1 -Norm Approximation nach dem 
Algorithmus von Barrodale und Roberts. Es werden, wie in 
Kapitel 9.1 besprochen, die Orienti er ungsun bekann t en in die 
Basis gebracht. Dann wird der normale Algori t hmus C vg I . 
Kapitel 6, Formeln C6.1)-C6.4)) wei t er verwende t . Die Matrix G 
wird in "random packing" abgespeichert, weit er s wird die 
LU-Zer I egung Cnach Kapitel 7) verwendet. Die Koeffizienten 
befinden sich immer blockweise im Kernspeicher und die Bioecke 
werden je nach Bedarf vom Disc-File nachgeladen. Eine Ze i I e 
besteht aus 16 Elementen und enthaelt alle Daten einer 
Beobachtung. In ihr stehen Informationen ueber die 
Beobachtung, Standpunkts- und Z i e I punk t snummer y Gewicht, 
gemessene Beobachtung, gemessene minus gerechnete Beobachtung 
Crechte Seite), Spaltenummern und Werte der Koeffizienten. 

Durch diese Abspe i cherungsar ten ergibt sich ein geringer 
Spe i cherpl atzbedarf . Fuer ein Netz mit 2500 Beobachtungen und 
200 Neupunkten Cpro Punkt ein gemessener Richtungssatz) werden 
30 K-Worte benoet igt . Das Programm benoetigt ca. 5 K— Worte. 

Die Rechenzeiten sind relativ hoch. Zum Beispiel braucht das 
Programm fuer ein Lagenetz mit 253 Beobachtungen und 78 
Unbekannten C15 Neupunkte, 48 Orientierungen) ca. 2 Minuten an 
CPU-Zeit CUNIVAC 1100/08). Davon entfallen auf Stufe 2 des 
Algorithmus von Barrodale und Roberts ca. 65% an CPU — Ze it. In 
Stufe 1 erreicht man schon eine Reduzierung der Absolutsumme 
der Verbesserungen um ca . 90% . 
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10. NUMERISCHE UNTERSUCHUNGEN FUER GEODAETISCHE LAGENETZE 



Fuer die numerischen Untersuchungen wurde ein Testnetz 
herangezogen. Es beinhaltet 253 Beobachtungen C 1 79 Richtungs-, 
74 Streckenbeobachtungen), 8 Festpunkte und 15 Neupunkte. Der 
mittlere Beobachtungsfehler fuer die Richtungen ist mit 6 
Neusekunden, fuer die Strecken mit 4 cm angegeben. Figur 10.1 
stellt eine Skizze dieses Netzes dar. Es wurden dam i t eine 
Reihe numerischer Untersuchungen durchgefuehr t / die im 
einzelnen besprochen werden. 

10. 1 Suche nach groben Fehlern (Ausreissern) 



Durch die veral Igemeinerten Medianeigenschaften ist die 
Ausgleichung nach der 1 -Norm sehr s t ab i I gegenueber groben 
Fehlern. In dem Testnetz wurden folgende grobe Fehler an den 
Beobachtungen und Neupunkten angebracht. 

a.) Richtung 85133-13134 im 1 . Buende I : 

- 1 Neum i nute 
b. D Standpunkt 26132: 

Z i e I punk tsver wechs I ung zwischen den Punkten 

85 133 und 9 133 
in allen drei Richtungsbuendeln 
c.3 Richtungen von 84161—26132 in allen drei Buendeln; 

I Neugrad 
d . 2) Streckenfehler: 

13134-9133 10.0 cm 
207164-5131 16. 1 cm 
60163-30164 10.0 m 
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TESTNETZ 
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A. 
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;- 38182 



89132 . : : ' ;; -' .*•' 

, : Q /:£•*.'. 85133 



'••'A 13134 
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26132- 



9133 \ 
ö;' 



5,162 '.)''&':•-• 
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;/A 207164 



£; 



6oi63 ;■<*■■■ 



45168 / .;•.'•*' 25161 ', .•'';. ■■■■■"" 42182 .' P - 



A 



•ov-; . ;.. . 



10 



38164 



Vi. 1162 :: '•■:.,;■.. /..V3I63 ,.;. 
•Qv-'-V ; ; :. ; .V/.- : i> • • • ' v 



\ 72187 p.' ; '■; 



"• :'•&' ;ö 65189 

: 55I8B 



A. 

46187 



10 KM 
f 1 



18183 



Flg. 10. ) 
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1 162-72187 10.0 cm 
q . ) Or i ent i erungsf eh I er von 1 Neugrad im 1 . Buende I von 

Standpunkt 65189 
f.) Fehler In den Naeher ungswer ten : 

9133 50 Meter in X-Richtung 

1 162 10 cm in X- und Y-Richtung 

Tabe I I en 10. 1 bis 10. 12 stel I en die Compu t er ausdrucke des 

Endergebnisses dar. Die 3 Sterne zeigen jene Beobachtungen an, 

die exakt den Wert nu I I bekommen CN i ch t bas i sver besser ungen) . 

Sie gehoeren demnach zu jenen Beobachtungen, die das optimale 
Netz aufbauen. 

Wie man aus den Tabe I I en ersehen kann, werden alle 
e i ngef uehr t en Fehler erkannt Cbuch gr oessenor dnungsmaess i g) . 
Vergleicht man die Ausgleichung mit den Ausreissern, mit jener 
ohne Ausreisser, so zeigt s i ch, dass zwischen den beiden 
Faellen kaum Unterschiede \ n den Verbesserungen Cnatuerli ch be i 
jenen, an denen kein grober Fehler angebracht wurde) bestehen. 
Bei den Richtungsverbesserungen differieren etwa 5X in der 
Gr oessenor dnung von ca. A Neusekunden. Die restl i chen bleiben 
gleich oder haben nur Unterschiede von ungefaehr 0. 1 
Neusekunden. Bei den Streckenverbesserungen differieren 
ebenf a I I s nur ganz wen i ge um etwa A cm . Die rest I i chen b I e i ben 
gleich gross oder weisen Aenderungen von einigen Mi I I imetern 
auf. Ebenso differieren die Zusch I aege bei den Neupunkten 
(wiederum jene, die keine groben Fehler haben) in beiden 
Fae I I en kaum. Von den 15 Neupunkten erha I ten 8 genau die 
gleichen Korrekturen und \ r\ den restli chen A Punk ten 

differieren die Korrekturen nur im Mi I I imeterbereich. 



123 



Die grossen Verbesserungen in 
Richtung 26161-84161 
Richtung 26161-46187 



R i ch tung 
S trecke 
St recke 
S trecke 
S trecke 



38102-10102 
26161-46187 
45160-84161 
51 162-65189 
207164-10102 



C 1 . Buende I ) 
C3.Buende I ) 
C 1 . Buende I ) 



treten auch schon in der Ausgleichung ohne die e i ngef uehr t en 
Ausre i sser in derselben Gr oessenor dnung auf und weisen somit 
auf schlechte Messungen in den eigenti i chen Beobachtungsdaten 
hin. 

Nur in Strecke 26132-9133 Cv=-0.098 m) hat die 
Verbesserung einen unterschied! i chen Wert im Gegensatz zur 
Ausgleichung ohne Ausre isser (v=— 0.038 rrO . Man muss aber 
bedenken, dass der Punkt 9 133 in X um 50.0 m falsch ist und 
somit die L i near i s i erungss t e I le relativ weit entfernt ist. 



Man erkennt hier die Stabi I i taet des Netzes, dass trotz 
grober Messfehler der restliche Netzteil unbee i nf I uss t bleibt. 



** ** + + ** MMMM** ****** ** ** ** ** ** ** * * **** 

♦LI - APPROXIMATION EINES EBENEN NETZES * 

** ************************** ** ************ 
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DIMENSIONEN : RICHTUNGEN IN NEUG RA D r STRECKEN IN KETER 
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10.2 Gew i chte 



Eine Aenderung des Verhae I Inisses der Richtungsgewichte zu 
den St reckengew i ch ton von 1 : 1 auf 1 :3 und 1 :9 bringt keine 
signifikante Aenderung im Ergebnis. 

Vorsicht ist al lerdings geboten, wenn man einzelne 
Messungen mit viel zu grossen oder viel zu klei nen Gew i chten 
e i nf uehr t . Ein Beispiel wurde gerechnet mit einer grob 
falschen Streckenbeobachtung mit sehr grossem Gewicht C 1 000) . 
Im Ergebnis hatte die grob falsche Strecke eine Verbesserung 
von Nu I I Cund war Nichlbas isvar i able) . Sinnvol I erscheint eine 
Gewichtswahl, die sich nach dem Perpendikel richtet. 



10.3 Fehler in der absoluten Orientierung des Netzes 



Bei freier Ausgleichung werden Fehler in den Parametern 
fuer die absolute Orientierung des Netzes (Verschi ebung, 
Masstab, Verdrehung) eindeutig erkannt. Es aeusser t sich dies 
durch systematische Zusch I aege in den Neupunkten und 
Orientierungen. Es bleiben aber die Verbesserungen und 
Korrekturen zu den Neupunkten und Orientierungen Cgemeint sind 
die Korrekturen mit Abzug der gemeinsamen Antei le von 
Verdrehung; Masstab und Verschiebung!) dabei ganz gleich wie bei 
der Ausgleichung ohne irdendwelche Fehler in den Parametern der 
absoluten Orientierung. Dies ist auch logisch, da ja die 
Absolutsumme der Verbesserungen dasselbe Minimum hat, wenn das 
ganze Netz einer Verdrehung, Verschiebung und Mass t absaender ung 
unterworfen wird. 
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Bei Zwangsausgleichung werden nur mehr Verdrehung und 
Verschiebung aufgedockt Cwenn z.B. al le Neupunkte eine 
Verschiebung gegenueber den Festpunkten aufweisen). Die Fehler 
werden wieder durch systematische Korrekturen in den Neupunkten 
ers i cht I i ch gemacht . Ein Masstabsfeh I er des ganzen Netzes 
Cwenn der Masstab der Streckenmessungen nicht mi t dem Masstab 
der Festpunkte vertraegi ich ist) bewirkt bei Zwangsausgleichung 
jedoch nur lokalen Einf luss. Sind zum Beispiel 2 Festpunkte 
vorgegeben, und ist der Masstab des ganzen Netzes zu gross, so 
werden nur bestimmte Neupunkte Cim Bereich 2 — 3 Punkte von den 
Festpunkten entfernt) zu den Festpunkten gezogen. Ebenso 
werden nur die Beobachtungen in diesem Te i I beeinflusst. Dem 
Netz wird eine lokale Verzerrung Cin der Naehe der Festpunkte) 
aufgezwungen. Neupunkte, die weiter von den Festpunkten 
entfernt sind, bleiben in dem Masstab, der durch die 
Streckenbeobachtungen vorgegeben ist, erhalten. 



10.4 Inhomogeni taeten in einem Netztei I 

Wieder werden Verschiebungen und Verdrehungen in einem 
Netztei I in den Zusch I aegen zu diesen Neupunkten erkannt. 



Tritt eine Masstabsinhomogenitaet in ei nem Netztei I auf, 
so bleibt die Beeinflussung durch diesen Fehler auch wieder nur 
lokal beschraenkt . Der Netzteil schrumpft oder dehnt sich Cje 
nach Masstabsfehler). Die Nachbar punk t e zum inhomogenen 
Netztei I werden nur ganz unwesentl ich beeinflusst, der uebr i ge 
Netztei I bleibt ganz gleich Cim durchgerechneten Beispiel hatte 
der uebr i ge Netztei I dieselben Verbesserungen und Zusch I aege 
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wie bei der Ausgleichung ohne Mass t abs i nhomogen i t ae O . Diese 
'endenz bleibt dieselbe, egal ob freie Ausgleichung oder 
Zwangsausgleichung. Bei Zwangsausgleichung kann jedoch der 
nhomogene Netzte i \, wenn er sich in der Naehe eines 

-esipunktes bef indet„ zu diesem noch zusaetzl ich hingezogen 
-er den . 

Die gleichen Merkmale zeigen s i ch, wenn zum Beispiel in 
einem Standpunkt al le gemessenen Strecken einen Masstabsfehler 
aufweisen. Hier ist jedoch eine starke Tendenz zu bemerken, 
■z e Strecken richtig zu verbessern. Wie auch zu erwarten war , 
Dekommt eine Strecke die Verbesserung nu I I CN i ch t bas i s) . Und 
zwar ist das dann eine jener Beobachtungen, die zur Festlegung 
aes Standpunktes notwendig war. Im folgenden Kapitel wird eine 
~aehere Erk I aerung dafuer gegeben. 



10.5 Auswahl der Beobachtungen fuer den optimalen Netzaufbau 

Die Berechnungen haben gezeigt, dass die in Kapitel 6 
Geschriebenen Eigenschaften des optimalen Netzaufbaues Ceine 
optimale Loesung steht im Zusammenhang mit einer Auswahl von 
Beobachtungen, die ausreicht das Netz widerspruchsfrei zu 
bestimmen!) viel exakter angegeben werden kann. Diese globale 
Aussage fuer das ganze Netz kann nun auf jeden einzelnen Punkt 
beschraenkt werden. 



Fuer jeden Punkt haben so viele Beobachtungen die 
Verbesserung nu I I Csind N ichtbas isvar i ab I e), wie zur 
widerspruchsfreien Festlegung dieses Punktes notwendig sind. 
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Es f uehr t auch folgende Ueber I egung zu diesem Ergebniss , Es 
sei das Optimum fuer alle Punkte des Netzes bis auf einen 
einzigen vorgegeben. So I I nun fuer diesen einen Punkt die 
Absolutsumme der Verbesserungen zu einem Minimum gemacht 
werden, dann muessen auch hier diese Eigenschaften des 
-optimalen Netzaufbaues gelten; es muessen genau so viele 
Beobachtungen die Verbesserung nu I I haben, wie zur 
einwandfreien Bestimmung des Punktes notwendig s i nd . D i ese 
Ueberlegung laesst sich natuerlich auch mit allen anderen 
Punkten des Netzes anstel I en . 

10.6 Beschr aenkung der Ausgleichung auf kleinere Te i Inetze 

Wie aus den bisherigen Untersuchungen folgte., bleibt der 
E i nf I uss von Fehlern lokal beschraenkt. Auch in der Praxis 
treten Ausreisser meist nur lokal in verschiedenen Gebieten des 
Netzes auf. Numerische Untersuchungen von kleineren Bereichen 
haben gezeigt, dass die Fehlersuche auch auf kleinere Te i Inetze 
beschraenkt bleiben kann. AI lerdings ist auf eine genuegende 
Ueberbest i mmung jedes einzelnen Punktes zu achten. Die 
Teilnetze sollen so aufgesplittert werden, dass die i nneren 
Punkte moegl ichst wenige Beobachtungen in andere Tei Inetze 
aufweisen. Vorsicht ist in den Randpunkten geboten, hier 
koennen durch die fehlenden Beobachtungen C die sich im 
angrenzenden Tei Inetz befinden) falsche Seh I uesse gezogen 
werden Cdies ist dann der Fa I I, wenn der Randpunkt eine zur 
Ausr e i sser I oka I i s i erung ungenuegende Ueberbes t i mmung hat). 



Guenstig ist es, wenn sich die Teilnetze im Randbereich 
mit ca. 2 Punkten ueber I appen . Grobe Fehler koennen aber in 
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Tel Inetzen gleich gut erkannt werden wie im ganzen Netz. 



10.7 Berechnung von Naeherungswer ten 



Sind bei der Ausgleichung nach der 1 -Norm die 
Naeherungswer t o al I er Neupunkte nur ganz grob bekannt Cetwa aus 
einer Karte mit einer Genauigkeit von ca. 1000 rrO, so ist eine 
Suche nach groben Fehlern nicht moeg I ich. Die Naeher ungss t e I le 
ist dann zu weit entfernt. Im speziellen gerechne ten Be i sp i e I 
CNaeher ungswer t e auf etwa 1000 m genau bekannt) wurden durch 
die Zusch I aege die Neupunkte nur auf etwa 20 m genau bestimmt. 
Auch in den Verbesserungen treten dann grosse Werte auf. Eine 
neue Iteration ergibt dann schon gute Wer te, in der auch grobe 
Fehler aufgedeckt werden koennen . 



10.8 Fehlersuche nach Ende von Stufe 1 

In den meisten Fae I I en haette zur Suche nach groben 
Fehlern schon Stufe 1 genuegt. Ungenuegende Ergebnisse sind 
nur bei Netzen, die Fehler in der absoluten Orientierung 
hatten, aufgetreten. Trotzdem ist der Abbruch nach Stufe 1 
nicht zu raten, da es i n E i nze I f ae I I en zu falschen Schluessen 
fuehren kann. 



10.9 Statistische Untersuchungen 



Um eine Aussage ueber das statistische Verhalten der 
Verbesserungen und Unbekannten zu bekommen, wurde das Monte 
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Carlo Verfahren angewandt. Die Beobachtungen wurden mit einem 
Zuf a I I sgener ator verfaelscht. Dio Zufall szah len wurden 
normalverteilt gewaehl l, fuer Ri chtungen mit e i nem 
Beobachtungsfehler von 4 Neusekunden und fuer die Strecken mi t 
einem mittleren Beobachtungsfehler von 4 cm CR i ch t ungsgew i ch t : 
Streckengewicht 1 : 3). Fuer jede Beobachtung wurden 100 
Zuf a I I szah I en generiert und das Netz nach der 1— Norm berechnet. 

Die Ausgleichung nach der 2-Norm (vermi t te I nd, mit 
Beobachtungsfehlern vor der Ausgleichung: Richtungen 6 
Neusekunden,, Stecken 2 cm bis 4 cm) ergab folgende Werte nach 
der Ausgleichung: 

mittlerer R i ch tungsf oh I er 4.2 Neusekunden 

mittlerer Streckenfehler 1.4 cm 

durchschnitt! i qher 

mittlerer Koordinatenfehler 1.6 cm 



Die Ausgleichung nach der 1— Norm ergibt fuer die 
Verbesserungen CBas i sver besser ungen ) und fuer die Zusch I aege zu 
den Unbekannten C Koord i na ten i nkremen te , Orientierungen) wieder 
eine empirische Norma I ver te i I ung . Figuren 10.2 bis 10.5 geben 
fuer 4 Beobachtungen die relativen Haeuf i gke i ten nach der 
1— Norm Ausgleichung an. Und zwar werden hier die relativen 

Haeuf i gke i len ohne die exakten Nu I Iwerte (Nichtbasis) 
aufgetragen Cbedingte Verteilung). Die relative Haeufigkeit 
der exakten Nu I I werte, die auftraten, sind extra als dicke 
Linie auf der Ordinate eingezeichnet. Es sind Richtungen und 
Strecken zwischen Neupunkten, bzw. zwischen einem Festpunkt und 
Neupunkt dargestel It. Die restl i chen Beobachtungen weisen 
keine signifikanten Unterschiede auf. Die Standardabweichungen 
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der Richtungen bewegen sich zwischen 3 und 5 C dur chschn I t t I Ich 
4.3) Neusekunden, die der Strecken zwischen 3.5 und 5 cm 
C dur chschn ittlich 4. 5cm). 



D i e Wahr sehe inlichkeit Coder exakter, die empirische 
relative Haeuf i gke i t ) , dass eine Verbesserung 
Nichtbasisvariable wird, ist unterschied! ich. Hier muss 
zwischen Beobachtungen, an denen nur Neupunkte betei I igt sind 
und Beobachtungen, an denen ein Neupunkt und Festpunkt 
betei I igt sind, unterschieden werden. Es ergeben sich folgende 
empirische Werte fuer die Wahrscheinlichkeit: 

Richtung Neupunkt - Neupunkt 20X-40X 

Richtung Neupunkt — Festpunkt 1 7%-37% 

Strecke Neupunkt - Neupunkt 46%-72% 

Strecke Neupunkt - Festpunkt 28%-55% 
Bei Beobachtungen zwischen Festpunkten und Neupunkten ist die 
Wahrscheinlichkeit etwas niedriger Cdies ist auch 
verstaendl ich, da ja der Festpunkt einen zusaetzl i chen Zwang 
auf die Beobachtung ausuebl) . Interessant ist die Tatsache, 
dass hier das hoehere Gewicht der Strecken zum Tragen kommt, im 
Gegensatz zu den Standardabweichungen, wo sich die Gewichte 
nicht signifikant auswirken. Die Standardabweichungen in den 
Koordinaten der Neupunkte sind etwas hoeher als beim Ausgleich 
nach der 2-Norm Cim Durchschnitt etwa 3cm). Dies war auch zu 
erwarten, denn bei direkten Beobachtungen mit Normalvertei I ung 
ist die Standardabweichung des Med i ans 1 .25 mal so gross wie 
die des Mittelwertes. 
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Richtung (Neupunkt —Neupunkt) 
(26132 - 9133) 



rel ati ve 
Häufigkeit 



1-Norm Verteilung 

Normal Verteilung 



507.-- 



Standardabweichung: 4. 1 
Mittelwert: -Q.6 CC 



cc 



-exakt null (Nichtbasis) 



-11.3 -6.7 

Min. Wert 



-6.Ü -3.3 -0.7 



Fig. 10.2 



Richtung (Neupunkt— Festpunkt) 
(84161 - 51162) 



2.0 4.6 7.3 

Max. Wert 



relative | 
Häufigkeit I 

1-Norm Verteilung 50v- 



Normal vertei 1 ung 



— 



-11.0 -7.7 

Min. Wert 



-4.5 



-1.2 



Standardabweichung: 4.5 
Mittelwert: 0.3°° 



cc 



; exakt null (Nichtbasis) 



V 



2.1 



5.4 8.6 11.9 

Max. Wert 



Fig. 10.3 
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Strecke (Neupunkt —Neupunkt) 
(1162 - 3163) 



relative 
Häufigkei t 



1-Norm Vertei 1 ung 



Normalvertei lung 



Standardabweichung: 4.9 cm 
Mittelwert: -0.4 cm 



5 0°/. 




exakt null (Nichtbasis) 



-12.8 -8.6 

Min. Wert 



8.3 12.6 

Max. Wert 



Fig. 10.4 



Strecke (Neupunkt —Festpunkt) 
(51162 - 26132) 



relative 
Häufigkei t 



1-Norm Verteilung 

Normal vertei lung 



5 0Vc 



. -exakt null (Nichtbasis 



r 



-15.0 -10.6 

Min. Wert 



T" 



Standardabweichung: 4.9 cm 
Mittelwert: -0.12 cm 



6.3 



- 1.9 



IX 



ü 



6.8 



I 

11.1 



Max. Wert 



Fig. 10.5 
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I . NUMERISCHE UNTERSUCHUNGEN ZUR PUNKTTRANSFORMATION 

UNTER MINIMIERUNG DER ABSOLUTSUMME DER RESTKLAFFUNGEN 



An Hand eines Punkthaufens von 12 Punkten wurden folgende 
Eigenschaften getestet: 

a.) Aufdeckung von groben Fehlern 

b . ) \^qt a I I gerne i nerte Med i ar\Q i genschaften 

c.) Konfigural ionen, bei denen ein Punkt die Restklaffung 

nu I I erhalten muss 
d.) Inhomogen i t ae t en eines Te i I es der Punkte (Masstab, 

Versch i ebung) 

ad a . ) Der in Tabelle 11.1 dargestellte Punkthaufen Qj 
Ci = 1 , . . . I 2 !) wurde einer Helmerttr ansf or ma tion auf den 

Punkthaufen P. unterworfen. 

i 

Vorher wurden in den Punkten Q. „ Q ? , Q. Q folgende grobe 
Fehler eingefuehrt: 

Q, : Klaffung : .846 m 

Richtung: 298.420 Neugrad 



Q ? : Klaffung: 1.151 



'2 



m 



Richtung: 340.580 Neugrad 



Q in Klaffung: 0.253 



m 



Richtung: 400.000 Neugrad 
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Pktnr . 


P 






Q 




1 


Y 


X 


Y 


X 


1000 


000 


9500 


000 


1000 


000 


9500 


000 


2 


10000 


000 


10500 


000 


10000 


000 


10500 


000 


3 


9500 


000 


6000 


000 


9500 


571 


5999 


354 


4 


9000 


000 


2500 


000 


9000 


296 


2499 


393 


5 


6000 


000 


500 


000 


6000 


139 


499 


628 


6 


2500 


000 


2500 


000 


2500 


296 


2499 


903 


7 


3000 


000 


8500 


000 


3000 


768 


8499 


864 


8 


6000 


000 


9000 


000 


6000 


807 


8999 


629 


9 


8000 


000 


5500 


000 


8000 


532 


5499 


472 


10 


4500 


000 


4500 


000 


4500 


453 


4500 


000 


1 1 


3500 


000 


5500 


000 


3500 


532 


5499 


825 


12 


5727 


273 


5863 


636 


5727 


.833 


5863 


287 



Tobe I ie 11.1 

Die Koordinaten der Punkte Q wurden so gewaeh I t , dass eine 
differentiöl Ie Drehung von 50 Neusekunden und eine Verschiebung 
in X— und Y— Richtung von 0.1 m gegenueber den Punkten P 
herrscht, ferner so angenommen, als ob die 
NaehGrungstransformat ion mit den grob falschen Punkten 1 und 2 
dur chgef uehr t worden waere. 



Wie aus dem Computerausdruck in Tabe I Ie 11.2 zu ersehen 
ist, stimmen die Transformationsparameter sehr gut mit den 
vorher s i mu I ierten ueberein. Ebenso die Groesse und Richtung 
der e i ngef uehr t en groben Fehler wird eindeutig aufgedeckt. 
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Fig. 11.1 
Skizze der in Tobe I le 11.1 angegebenen Punkte 

M 1 = 100000 



PN R 



OY OX R VY VX RhO DR TAU LAfÖDA PI 



1 .CirtOU .OiJOCü .64543 -.*45:6 -.02149 . 8 4 5 4 3 .IHUlCU 298.38154 298.38154 1 .3 03 00 ♦ 

2 .00000 .0 00 CO 1.15032 -.92436 .66470 1.15032 .O'VjOh 340.5H699 3 40. 58699 l.DUOOO * 

3 .57100 -.6«6Chi ,üü!)?6 -.liCOie -.00018 .Ofkl26 .OCUCfl 251.066(13 251.06611 1.3 00 03 ♦ 

4 .29600 -.6070.1 .00050 -.OOOU -.UC013 .00050 .POuOO 282.972C5 262.97199 1 .3 03 03 

5 .I39C0 -.372,00 .0004-9 -.9C029 -.''0039 . :. CO 4 9 .OPiiGO 240.8J983 24J.H3976 .7470J ♦ 
m m \ m m ? 40 .8399 1 .2 53 Oü • 

g- 6 .29600 -.097UO . JOJ22 .ü^OOl -.U0a2I .UG022 .'JOOGO 188.21422 1 tri .2 1 4 3 3 .11929 » 

ß ....*. . • 1 £B .2 1 40 5 .fth.l 71 

^ 7 .768CO -.13600 .OT124 .00116 -.00046 .1)0124 .00003 124. 21525 124 .21530 1 .3 U3 CO ♦ 

8 .8C7CO -.3/100 .JJM68 .00063 -.00002 .00-168 .00000 107.46380 107.48387 .37342 • ^ 

I 07 . 4 R 3 7 2 .92958 • 



l\) 



9 .532CO -.52800 .r»OÖ23 .00017 .0 0:1 1 5 . -if.-i 23 .OOOCO 54.4817M 54 .48162 .81?<U • 
5 4.4 8 194 -l £7 5 9 • 

10 .45300 .OOL'CG .25357 -.00007 .2535/ .25357 .00000 399.98255 399.98255 1.3U300 

11 .532CO -.17500 .'J^054 .OGC54 .00002 .<<G:i54 .00000 97.94140 97.^4150 1,3 03 00 ♦ 

12 .56027 -.34936 .000^2 .10009 .00041 .00042 .00000 14.44352 14.44358 1 . 3 üü CO • 

TR AOSf OP f'A TI 0\'9 AR AM El £»: 

MA ^S T AD = - . 00 00 80 7 

VERDREHUNG - -.07 84 74 1 [BOGENMASS • 1000 ) 

VERSCHIEBUNG IN X - -.0991981 [m] 

VERSCHIEBUNG IN Y ^ -. 099572 7 [m] 



PUR 



ü Y 



D.< 



V Y 



V > 



RHC 



CR 



TAU 



L AtfßDA 



PI 



Q 

CT 

Q 



U) 



! .00000 .0 03 00 .84543 - . «* i Z <j — • ' i P 1 «4 V . b 4 5 4 3 ,n OOOO 2 98 . 3 b 1 54 2 93 . 3 8 3 5 4 1.3C3 0Ü ♦ 

2 ~56.2<J9<n 41.*'> r >37 71.15032 -57.17 37 / 42.35O07 75.15032 .00000 340.58699 343.58699 1 .3 03 CO ♦ 

3 .571 OC -.64600 .0 00 2* -.00016 -.'J'ÜIIB . ') Oü 2 6 .flCJiJ't 2M.U66IM 251.06611 1.3 03 00 * 

4 .2960C -.6i}7 r )0 .'JCiJä't -.OCL^e - . i ' » ; f t i 3 .00050 .QOOÜO 262.97205 2F.?.v71W 1 ,3 U3 OJ 

5 .13900 -.37200 .:k\'*49 -.'.-0029 -.0 0;» 39 .1/0049 ..TTiC'"» 2 4 n . H 3 9 P 3 240. "3976 . 7j9 95 • 
, . . . . 2 ** 3 . h 3 9 9 l .26305 ♦ 

6 .2960-1 -.097CO - 1 1 i : 2 2 .^0004 -,0«3'l2L .00022 .00^00 1^8*22422 1 PK . ? UO r ; .bßI35 » 
I 88 .2 1 «3 S -l 1 B 65 

/ .76800 -,13600 .:>-31P«4 .MIU -.U0H46 .OCi:«! .nCHOO 124.21525 124.21530 2.30303 • 

& .60700 -.37100 .0006« .00068 -.00008 .0006d .30000 107.48379 107.48387 .3 96 90 * 
. . . . . . 107.4R372 .90310 * 

9 .53200 -.528CO .30023 .IMÜ17 .00015 -IV)M2 3 .0 00 00 54.48178 54.48162 .81545 • 
. . 54.48194 . 1 d 4 5 5 * 

10 .'55300 .On:i."!0 .253^7 -.n0.'!')7 .25357 .25357 .00000 399.9ä255 399.98255 1 .3 C3 03 ♦ 

13 .53201 ~.17 r >0<i ,00054 .00054 .00002 .00054 .00000 97.94140 97.94150 1 .3 03 00 • 

12 .56027 -.3*4936 .00')m2 .'IUüÜ9 .'JOH41 .0-1042 .OOOQO 14.44352 14-4435Ü 1 .3 03 CO * 



CD 



TK Mi 5F CR MA T I OM SP AR AM E T £ R i 



MA SS TAB 

VL RiREHUr*G 

\Jt PS CM IE BUN'C I N X 

VE RS CM 1F BUNC I N Y 



-.00 00 80 7 

-.07 64 74 1 [BOGENMASS * 1 OOOJ 

-.09 9t 98 1 [M] 

-.0995727 [M] 
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Die Iterationen bei der Berechnung wurden solange 
dur chgef uehr t , bis TAU und LAMBDA im Sekundenbereich 
uebereinstimmen. R , V Y , VX in Tabelle 11.2 sind die aus den 
Transformationsparametern berechneten Restklaffungen und 
Verbesserungen. RHO ist die Res tk I af f ung, die sich aus der 
Loesung des I inearen Programmes ergibt. DY und DX sind die 
Koordinatendifferenzen vor der Transformation. In den Punkten 
5, 6, 8, 9 wird der Kreis in der Richtung TAU durch 2 gleiche 
Tangenten angenaeher t . Die Summe der PI CID ergibt aber wieder 
1 . 

ad b.) Im Punkt 2 wurde die urspruengl i che K I af f ung in 
Richtung von TAU um 70.0 m erweitert. Hier zeigt sich die in 
Kapitel 8.7 besprochene verallgemeinerte Medi ane i genschaft . 
Der dazugehoer i ge Compu t er ausdruck ist in Tabe I le 11.3 
dargestellt. Es zeigt sich, dass saemtliche Restklaffungen und 
Transformationsparameter gleich wie bei Fa I I a.) bleiben. Nur 
die Restklaffung im Punkt 2 steigt um genau 70.0 m. Auch die 
Richtungen der Restklaffungen bleiben erhalten. 

ad c.) Zur numerischen Testung einer Konf i gur a t i on, in der 
e i n Punkt d i e Restk I af fung nu I I erhae I t, wurde der Punkt 1 2 
weit hinausgeschoben CY = 100000.000, X = 100000.000). Es 
zeigt sich, dass dieser Punkt tatsaechl i ch, nach den in Kapitel 
8.10 besprochenen Eigenschaften, die Restkl äff ung nu I I er hae I t . 

ad d.) Inhomogenitaeten eines Teiles der Punkte (z.B. im 
y asstab oder eine Verschiebung) werden in den Restklaffungen 
eindeutig erkannt. 



15 



12. ZUSAMMENFASSUNG 



Die Ausgleichung durch Minimierung der Absolutsumme der 
Verbesserungen ( I -Norm) eignet sich vorzuegl ich zur Aufdeckung 
grober Fehler in den Beobachtungen und Unbekannten. Die 
Loesung erfolgt auf dem Wege der I inearen Optimierung, Der 
Algorithmus von ßarroda I e und Roberts bietet hier das bisher 
schnei Iste und am wenigsten Speicherplatz benoetigende 
Verfahren an. Die Anwendung von Drei ecksz erlegungen bei 
Auf I oesung von den Gleichungsystemen bietet grosse numerische 
Stabil itaet. Das Abspeichern saemtl icher Vektoren und Matrizen 
in "random pack ing" bringt sehr grosse Ersparnisse im 
Speicherplatzbedarf, man muss jedoch groessere Rechenze i ten 
dam i t in Kauf nehmen . 

Die Standardabweichungen in den Beobachtungen sind etwas 
groesser als bei der Ausgleichung nach der 2— Norm und 
entsprechen in etwa den Beobachtungsfehlern vor der 

Ausgleichung. Eine wichtige Tatsache ist, dass die 
Beeinflussung durch grobe Fehler immer nur lokalen Charakter 
hat. Daher bietet sich die Moeglichkeit an, z.B. grosse Netze 
in kleinere, sich ueber läppende Teilnetze zu unterteilen und 
diese einzeln auszugleichen. 

Die Abso I u t summenm i n i m i er ung kann nicht nur zur 
Fehlersuche in geodaeti sehen Netzen, sondern auch auf alle 
anderen Gebiete der Geodaes i e angewandt werden. Die 
He I mer ttransformation wurde schon in dieser Arbeit bespr ochen . 
Sehr gute Ergebnisse erzielte auch Schmid, H. C1980) bei der 
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3uende I ausg I e i chung nach der l~Norm. Auch hier lassen sich 
w i oder grosse pho t ogr amme tr i sehe Bioecke in Te i I b i oecke 
zerlegen. Meiss I, P. C1980) hat mit Hi I f e der Graphentheorie 
ein Verfahren zur 1 -Norm Ausgleichung eines Hoehennetzes 
er stellt. Der Algorithmus ist eine Spezialisierung der Methode 
von Barrodale und Roberts. Additionen und Subtraktionen sind 
dabei die einzigen Rechenoperationen, die benoetigt werden. 



Die Ausgleichung nach der 1 —Norm so I I nicht als 
Ausgleichung im herkoemm I i chen Sinne verstanden werden, sie 
so I I und kann auch die Ausgleichung nach der 2 — Norm nicht 
ersetzen. Vielmehr stel It sie ein Verfahren zur Aufdeckung 
grober Fehler in den Beobachtungsdaten und Unbekannten dar . 
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